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The ratio of the proton elastic electromagnetic form factors, Gep/Gmp, was obtained by measuring 
Pt and Pe, the transverse and longitudinal recoil proton polarization components, respectively, for 
the elastic ep ep reaction in the four-momentum transfer squared range of 0.5 to 3.5 GeV'^. 
In the single-photon exchange approximation, the ratio Gep/Gmp is directly proportional to the 
ratio Pt/Pt- The simultaneous measurement of Pt and Pi in a polarimeter reduces systematic 
uncertainties. The results for the ratio Gep/Gmp show a systematic decrease with increasing , 
indicating for the first time a definite difference in the distribution of charge and magnetization in 
the proton. The data have been re-analyzed and systematic uncertainties have become significantly 
smaller than previously published results. 

PACS numbers: 25.30.Bf, 13.40.Gp, 24.85. +p 



I. INTRODUCTION 

One of the fundamental goals of nuclear physics is to 
understand the structure and behavior of strongly inter- 
acting matter in terms of its basic constituents, quarks 
and gluons. An important step towards this goal is the 
characterization of the internal structure of the nucleon; 
the four Sachs elastic electric and magnetic form factors 
of the proton and neutron, Gsp, Gmp, Geu and Gmu, 
are key ingredients of this characterization. The elastic 
electromagnetic form factors are directly related to the 
charge and current distributions inside the nucleon; these 
form factors are among the most basic observables of the 
nucleon. 

The first direct evidence that the proton has an inter- 
nal structure came from a measurement of its anomalous 
magnetic moment 70 years ago by O. Stern [J]; it is 2.79 
times larger than that of a Dirac particle of the same 
mass. The first measurement of the charge radius of the 
proton, by Hofstadter et al. 0, yielded a value of 0.8 fm, 
quite close to the modern value. 

The theory that describes the strong interaction be- 
tween quarks and gluons is Quantum Chromodynamics 
(QCD). Perturbative QCD (pQCD) makes rigorous pre- 
dictions when the four-momentum transfer squared, Q^, 
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is very large and the quarks become asymptotically free. 
It is not known precisely at what value of pQCD may 
start to dominate; however, expectations are that this 
will not occur until is at least in the tens of GeV^ . 
Predicting nucleon form factors in the non-perturbative 
regime, where soft scattering processes are dominant, is 
very difficult. As a consequence there are many phe- 
nomenological models which attempt to explain the data 
in this domain; precise measurements of the nucleon form 
factors are necessary to constrain and test these models. 
Only the magnetic form factor of the proton, Gmp, is 
known with very good accuracy in this region. The elec- 
tric form factor, Gsp, was not well measured beyond 
of 1 GeV^ before this experiment. Both Geu and Gmu, 
the electric and magnetic form factors of the neutron, re- 
spectively, were also poorly known at any value until 
recently. New measurements of Gmu at Jefferson Lab 
up to Q^—4:.8 GeV^ will bring the knowledge of this 
form factor to comparable levels of accuracy as for Gmp- 
For the neutron electric form factor, two new Jefferson 
Lab experiments 0, 0| have extended the range to 
1.5 GeV^, and two approved experiments fT'^'q will soon 
extend the range to 4.3 GeV^, with an accuracy com- 
parable to that of the three other form factors. 

The electromagnetic interaction provides a unique tool 
to investigate the internal structure of the nucleon. The 
measurement of electromagnetic form factors in elastic, 
inelastic, and structure functions in deep inelastic scat- 
tering of electrons and muons, has been a rich source of 
information on the structure of the nucleon. 

In the single virtual photon exchange approximation 
for elastic scattering, the hadron current operator can be 
expressed in terms of two form factors: -Fi, the Dirac 
form factor, and F2, the Pauli form factor. These form 
factors and the Sachs electric and magnetic form factors 
are related according to: 

Ge = Fi- TKF2 and Gm = Fi + kFj, (1) 

where r = Q^/AM^, k is the anomalous magnetic mo- 
ment and M the mass of the proton. In the limit 0, 
Gep = 1, Getl = 0, Gmp = Mp, and Gmu = fJ-n, where fj.p 
and fj,n are the nucleon magnetic moments. In the Breit 
frame, Ge and Gm are the Fourier transforms of the 



3 



charge and magnetization distributions in the nucleon, 
respectively. 



A. Previous Gsp Measurements Using the 
Rosenbluth Separation Method 

Both the elastic cross section and the polarization ob- 
servables of the elastic ep reaction can be expressed in 
terms of either the Sachs or the Dirac and Pauli form fac- 
tors. These form factors are Lorentz scalars and depend 
only upon Q^, the four- momentum transfer squared of 
the reaction. A complete separation of the electric and 
magnetic terms is evident in the cross section formula 
when the Sachs form factors are used. It is then possi- 
ble to obtain both and separately, using the 
Rosenbluth method |9l Il0||. In the one-photon exchange 
approximation, the cross section in terms of the Sachs 
form factors can be expressed as: 
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Ee cos ^ 
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where e = [l -I- 2(1 -I- r) tan^(^)] ^ is the longitudinal 
polarization of the virtual photon, with values between 
0< e <1, Eiieam and are the energies of the incident 
and scattered electron, respectively, and 9e is the electron 
scattering angle in the laboratory frame. 

Figure 1 shows previous results of Gsp and Gjv/p ob- 
tained by Rosenbluth separations, plotted as the ratios 
Gep/Gd and Gmp/ IJ-pGo versus Q^, up to 6 GeV^. 
Here Go = (1 + Q'^I'^d)~'^ is the dipole form fac- 
tor, with the constant m|, empirically determined to be 
0.71 GeV^ For < 1 GcV^ the uncertainties for both 
Gep and Gup are only a few percent, and one finds 
that Gupl [i-pGn — Gep/Gd — 1- For Gep above 
— 1 GeV'^, the large uncertainties and the scatter in re- 
sults between different experiments, as seen in Fig. 1, 
illustrate the difficulties in obtaining Gep by the Rosen- 
bluth separation method. In contrast, the uncertainties 
for Gmp obtained from cross section data with the as- 
sumption Gpn = GmpIp^p, remain small up to = 
31.2 GeV^T?!. In Eq. ^ the Gmp part of the cross sec- 
tion, which is about /^^ times larger than the Gep part, is 
also multiplied by r; therefore, as increases, the cross 
section becomes dominated by the Gup term, making 
the extraction of Gep more difficult by the Rosenbluth 
separation method. 
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FIG. 1: World data prior to 1998 for {&)Gep/Gd and (b) 
Gmp/ t^pGn versus Q^. Refs. Litt et al. ]0|A, Berger et al. 

n , Price et al. I1J». Bartel et al. ll4|o. Walker et al. 
[l5| *. Andivahis et al. (liqlO and Sill et al. p^*. 



are able to obtain the electric form factor Gep even when 
it is very small. 

For one-photon exchange, in the ep — > ep reaction, 
the scattering of longitudinally polarized electrons results 
in a transfer of polarization to the recoil proton with 
only two non-zero components, Pt perpendicular to, and 
Pi parallel to the proton momentum in the scattering 
plane. For 100 % longitudinally polarized electrons, the 
polarizations are [Tl Il9ll20ll2l| : 



loPn = 



(3) 



loPt = -2^Til + T)GEpGAiptaTij (4) 



loPe = J^{Ebearn+E,)yMl+^Gl,ptl,Il^^{5) 

where Iq is proportional to the unpolarized cross section 
and is given by: 



B. Polarization Transfer Method 

The proton form factor ratio Gep/Gmp can be ob- 
tained from polarization observables of the ep ep or 
ep — > ep reaction, the recoil proton polarization transfer 
coefficients or the beam-target polarization asymmetry, 
respectively. Both reactions contain an interference term 
proportional to GepGmp, hence polarization experiments 



-^0 — + -Gf^fp (6) 

Eqs. Q and Q show that loPt and IqPi are pro- 
portional to GepGmp and G\ip, respectively. Together 
these equations give: 

Gep _ Pt {Ebeam + E^) 9e 

Gmp Pi 2Mp 2 
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If only the polarization components Pt and are mea- 
sured, as was the case in this experiment, then from Eqs. 
Q and the form factors Gep and Gmp cannot be 
obtained separately, only their ratio can be determined. 
To obtain Gep and Gmp separately, /q in Eq. © must 
be obtained from cross section measurements. 

The ratio Gep/Gmp is obtained from a single mea- 
surement of the two recoil polarization components Pt 
and Pi in a polarimeter, whereas the Rosenbluth method 
requires at least two cross section measurements made at 
different energy and angle combinations at the same Q^. 

The recoil polarization method was first used in elec- 
tron scattering experiments to obtain the neutron form 
factors in the '^H{e,e'fi)p reaction and to measure 
the form factor ratio Gep/Gmp in ep ep for the free 
proton j23|,|2J|, as well as in the ^H{e, e'pjn reaction for 
the proton in the deuteron at small Q^-values |25l |. 

For completeness we mention here that a small nor- 
mal component P.^'^ is induced by two-photon exchange 
mechanism, independent of beam polarization. The ob- 
servables Pt and Pe of this experiment are entirely due to 
polarization transfer, and the analysis method described 
in this paper allows complete separation of helicity de- 
pendent and helicity independent polarization compo- 
nents. 

In this paper, we present the Gep/Gmp ratios, pri- 
mary results from this experiment, obtained at Jefferson 
Lab using the recoil polarization method described here. 
The experimental setup, in particular the focal plane po- 
larimeter (FPP), is described in part II. The data analysis 
is presented in part III; this part also includes a discus- 
sion of the FPP calibration, the secondary results of the 
experiment, which are independent measurements of an- 
alyzing powers at ten proton energies between 0.244 GeV 
and 1.795 GeV. Part IV includes the main results of the 
experiment: the ratios Gep/Gmp at 0.5 GeV^ < < 
3.5 GeV^, and an analysis of systematic uncertainties. A 
discussion of theoretical calculations as compared to the 
data is in part V and conclusions are presented in part 
VI. 



II. THE EXPERIMENT 

The combination of high energy, current, polarization, 
and duty factor, unique to the Continuous Electron Beam 
Accelerator Facility (CEBAF) of the Thomas Jefferson 
National Accelerator Facility (JLab), makes it possible 
to investigate the internal structure of the nucleon with 
higher precision than ever before. In this experiment, 
we have measured the polarization transferred to the re- 
coil proton, with a longitudinally polarized electron beam 
scattered by an unpolarized hydrogen target. 

The experiment was performed in Hall A at JLab. The 
longitudinal and transverse polarizations of the outgoing 
proton were measured for the ep ep reaction, in a 
range of from 0.5 GeV^ to 3.5 GeV^. The beam en- 
ergy ranged from 0.934 GeV to 4.091 GeV. For the five 



highest data points, a bulk GaAs photo-cathode ex- 
cited by circularly polarized laser light produced beams 
with polarization of ~0.39 and currents up to ^115 /iA; 
the sign of the beam helicity was changed at the rate of 
30 Hz. For the lower data points, a strained GaAs 
crystal was used and typical polarizations of ^^0.6 were 
achieved with currents between 5 fiA and 15 /iA; the sign 
of the beam helicity was changed at the rate of 1 Hz. 
The beam polarization was measured periodically with 
a Mott polarimeter in the injection line, and a M0ller 
polarimeter 26] in Hall A |27l |. 

The Hall A M0ller polarimeter uses magnetized fer- 
romagnetic supermendure foils as a polarized electron 
target. The scattered electrons are detected in co- 
incidence in the M0ller spectrometer in the range of 
75° < OcM < 105°. The M0ller spectrometer con- 
sists of three quadrupoles and a dipole magnet to bend 
scattered electrons toward the detector. The detec- 
tor contains two identical modules for coincidence mea- 
surements; each module consists of a plastic scintilla- 
tor and four blocks of lead glass. The M0ller scat- 
tering cross section depends on the beam and the 
M0ller target polarizations, Pe,i and Ptarg,i, respectively, 

CT CX (1 -I- I]i=X,Y,z(^" Ptarg,i Pe,i)), whcrC i= X,Y,Z 

defines the projections of polarization. The analyzing 
power An depends on the scattering angle Oqm and has 
its maximum at 9cm = 90°. Statistical uncertainty 
varies between 0.2 % and 0.8 % for each measurement. 
Total relative uncertainty of the beam polarization mea- 
surement is < 3 %, when systematic and statistical un- 
certainties are combined. In this experiment, the beam 
helicity cancels in the ratio Pt/Pi, so strictly speaking 
measurement of the beam polarization is not needed, but 
the beam polarization was measured periodically to en- 
sure that the beam was polarized and also to allow FPP 
calibration as explained in section III D. 

The beam current was monitored continuously during 
the experiment using resonant (RF) cavities. The beam 
current monitor (BCM) in Hall A consists of an Unser 
psj l monitor sandwiched between two RF cavities. The 
Unser monitor provides an absolute measurement of the 
current; the RF cavities are calibrated relative to the 
Unser monitor periodically. Both components are en- 
closed in a box to shield them from stray magnetic fields 
and for temperature stabilization. 

Beam position and direction at the target were deter- 
mined from two beam position monitors (BPM) located 
at a distance of 7.524 m and 1.286 m upstream of the 
target position during this experiment. Each BPM is a 
cavity with a four-wire antenna with wires positioned at 
±45° from the horizontal and vertical. The relative po- 
sition of the beam on the target can be determined to 
about 100 fim for currents above 1 fj,A by using the tech- 
nique of difference-over-sum between the signals from the 
four antenna wires. To obtain the absolute position of 
the beam, the BPMs are calibrated with respect to wire 
scanners which are located close to each of the BPMs at 
7.353 m and 1.122 m from the target. The wire scanners 
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are surveyed with respect to the hall coordinates. The 
beam position is recorded for every event. 

The cryogenic target contained three loops. Each loop 
included one 15 cm and one 4 cm aluminum cell; both 
cells have a diameter of 6.35 cm. The sidcwall thickness 
of each cell was 178 /im, and entrance and exit window 
thicknesses were 71 /im and 102 /im, respectively. As this 
experiment required only a liquid hydrogen target, only 
loop 3 was used; the other two loops were filled with he- 
lium gas at 0.12 MPa to save cooling power. The nominal 
temperature and pressure for the liquid hydrogen target 
during the experiment were 19 K and 0.17 MPa, respec- 
tively. The target density decreased by about 5 % [23 at 
an incident beam current of 120 /xA compared to its den- 
sity at 10 /iA (measured in an earlier experiment). The 
target assembly was housed inside a scattering chamber. 

The scattering chamber in Hall A is divided into three 
sections. The vacuum in all three sections is maintained 
at a level of 0.13 mPa. The bottom section is fixed to the 
Hall A pivot and the top part of the scattering chamber 
contains the target's cryogenic plumbing. The middle 
section of the chamber has an inner diameter of 103.7 
cm, a wall thickness of 5 cm of aluminum and a height of 
91 cm. The entrance and exit beam pipes are connected 
to this section. The scattered particles go through alu- 
minum exit windows 18 cm high and 406 /xm thick to the 
entrances of two high resolution spectrometers (HRSs). 

In order to reduce the heat deposition in a very small 
area of the target from an intense electron beam and 
to minimize corresponding target density changes, the 
beam was rastered before it strikes the target. The fast 
rastering system is located 23 m upstream of the target. 
The rastering system contains two sets of magnets, one 
to deflect the beam vertically, and the other to deflect 
it horizontally. The magnetic field varies sinusoidally at 
17.7 kHz in the vertical direction and 25.3 kHz in the 
horizontal direction. The typical rastered beam spot size 
at the target was « 3.5x3.5 mm^. 

A box located at the entrance window of each spec- 
trometer can contain three movable collimators. The up- 
per collimator is a stainless steel 5 mm thick sieve slit; 
it is used to study the optics of the spectrometers. The 
middle collimator is made of tungsten, and is 8 cm thick, 
6.29 cm wide, and 12.18 cm high, and is located at a dis- 
tance of 110.9 cm from the target. The bottom position 
is empty and performs no coUimation and is the one used 
in this experiment; the coUimation is then defined by the 
aperture of the magnetic elements of the HRS. The space 
between the exit window of the scattering chamber and 
the entrance window of each HRS consists of 20 cm of 
air. 

Elastic ep events were selected by detecting scattered 
electrons and the recoiling protons in coincidence, us- 
ing the two identical HRSs of Hall A. Each spectrometer 
consists of three quadrupoles and one dipole. The con- 
figuration is QQD„Q, two quadrupoles followed by an 
indexed dipole (n=-1.25), and a quadrupole. Both spec- 
trometers are designed for point-to-point focusing in the 



dispersive, and mixed focusing in the non-dispersive di- 
rection. The front quadrupole Ql has a magnetic length 
of 0.941 m and an inner radius of 0.15 m, and it is focus- 
ing in the dispersive direction. The quadrupoles Q2 and 
Q3 provide focusing in the non-dispersive direction and 
they both have an inner radius of 0.30 m and a magnetic 
length of 1.82 m. All three are superconducting, cos(20) 
quadrupoles with an outside cylindrical magnetic field 
return iron yoke. The quadrupole fields are monitored 
using Hall probes and Gauss- meters. 

The dipole in the HRS has a superconducting coil and 
warm iron configuration with shaped pole-faces and field 
gradient to help focus in the dispersive direction; its mag- 
netic field deflects the particles in a vertical plane by 45°. 
The bending radius of the dipole is 8.40 m with a central 
gap of 0.25 m and an effective length of 6.6 m. The nom- 
inal maximum central trajectory momentum is 4 GeV/c, 
the momentum resolution is of the order of 10~*, and the 
momentum acceptance is ± 5 %. The fleld in the dipole 
is monitored and regulated by an NMR probe; it is stable 
at the 10"^ level. 

The focal plane detector assembly for each spectrom- 
eter is enclosed in a metal and concrete shield house to 
reduce the background radiation. Both detector systems 
contain two vertical drift chambers (VDC), and scintil- 
lator arrays called SI and S2. In addition, the electron 
detector system contains a gas Cerenkov detector, and 
lead-glass arrays used as pre-shower and shower detec- 
tors; the hadron detector package contains an aerogel 
Cerenkov, a gas Cerenkov and the focal plane polarime- 
ter (FPP). The assembly of the hadron arm detectors is 
shown in Fig. [3 In this experiment only the VDCs and 
SI and S2 detectors were used on the electron side, and 
the VDCs, SI, S2 and the FPP on the hadron side. 

The two VDCs, installed close to the focal plane of each 
HRS, give precise reconstruction of positions and angles. 
The central ray of the spectrometer passes through the 
center of each VDC at 45° to the vertical. The two VDCs 
are separated by 33.5 cm. The active area of each one 
is 211.8x28.8 cm^ with two wire planes at 45° to the 
dispersive direction and perpendicular to each other. All 
VDCs are operated at a high voltage of 4 kV, and the 
gas mixture used in these chambers is 62 % argon and 
38 % ethane. The position resolution of each plane is « 

100 filR. 

In both HRS detector systems there are six scintillator 
paddles in plane SI and six in S2. Each paddle is seen 
by one photomultiplier at each end. The paddles in both 
planes are oriented such that they are perpendicular to 
the spectrometer central ray. The distance between the 
two planes is 1.933 m on the electron side and 1.854 m 
on the hadron side. The active area of the SI plane is 
170 X 36 cm^, and of the S2 plane 212x60 cm^. The 
thickness of each paddle in both planes is 0.5 cm. The 
scintillator paddles of both planes overlap by 0.5 cm to 
achieve complete coverage of the focal plane area. 

The trigger for both spectrometers is similar and is 
formed from a coincidence between the signal from two 
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chambers 



FIG. 2: Schematic of the hadron arm detector package in- 
cluding the polarimeter. 



scintillator planes, SI and S2. The first requirement is 
to form a coincidence between the left and right signals 
from each paddle in the SI and S2 planes. The time res- 
olution per plane is about 0.3 ns (la). These coincidence 
signals are fed into a memory lookup unit (MLU) which 
is programmed to form a second trigger, called "S-Ray" 
(trigger Tl for electron HRS and trigger T3 for hadron 
HRS), by requiring that the paddles that fired in the Sl- 
and S2 planes belong to an allowed hit pattern. The al- 
lowed hit pattern requires that if paddle |N| is fired in 
the SI plane, then in the scintillator plane S2 a signal 
must come from paddle |N-1| or |N| or |N-|-1| or the over- 
lap between two of those. Each spectrometer also has 
a "loose" trigger: for the electron arm it is formed by 
requiring that signals from two out of three detectors be 
present, Sl-plane, S2- plane and the Cerenkov detector; 
the hadron arm loose trigger requires signals from just 
Sl-plane and S2-plane. These loose triggers are used to 
obtain detector efficiencies. Finally, the S-ray trigger sig- 
nals Tl and T3 form the "coincidence" trigger T5 for the 
experiment. These five different triggers, Tl to T5, are 
sent to the trigger supervisor (TS) and are also counted 
by scalers. 

The TS was designed and built by the CEBAF online 
data acquisition (CODA) group. The functions of the TS 
include interface between the hardware trigger electron- 
ics and the computer data acquisition system, producing 
a computer busy signal that is used to calculate the com- 
puter dead-time, and pre-scaling of the trigger inputs Tl 
to T5. 

The data acquisition system was entirely developed at 
JLab by the CODA group 30] . In this experiment CODA 
was running on a single Hewlett-Packard (IIP9000) com- 
puter. The two main tasks of CODA are to transmit 
information from the detectors to the computer via read 
out controllers (ROC) and build events by collecting data 
from all the ROCs. The data are stored on a hard disk 



temporarily for on-line analysis to monitor the exper- 
iment and then transfered to tapes in the JLab mass 
storage system to be used later for final off-line analysis. 

A. Focal Plane Polarimeter 

Polarization experiments have become increasingly im- 
portant in the study of nuclear reactions. Focal plane 
polarimeters were standard equipment at intermedi- 
ate energy proton accelerators, such as LAMPF [sif . 
TRIUMf|32|, SATURNE m, and PSI |3|. Experi- 
ments at these facilities have demonstrated the sensitivity 
of spin observables to small amplitudes. Similar consid- 
erations have more recently led to the development of 
proton polarimeters for use at electron accelerators, such 
as the MIT-Bates laboratory [3^1 and at the Mainz Mi- 
crotron j36j. 

The FPP in Hall A at JLab was designed, built, in- 
stalled and calibrated by a collaboration of Rutgers Uni- 
versity, the College of William and Mary, Norfolk State 
University, the University of Georgia, and the University 
of Regina [13 ■ 

1. Physical Description of Focal Plane Polarimeter 

The FPP is a part of the hadron detector package of 
the Hall A high resolution spectrometer. As shown in 
Fig. 121 the polarimeter is installed downstream from the 
focal plane VDCs; it is oriented along the mean particle 
direction in the focal plane area, at 45° to the vertical. It 
consists of two front detectors to track incident protons, 
followed by a carbon analyzer and two rear detectors to 
track scattered particles. 

The four tracking detectors of the FPP are drift cham- 
bers made of straw tubes; the straw tube design is based 
on the one used for the EVA detector at Brookhaven 
National Laboratory The four drift chambers con- 
tain a total of 24 planes of straw tubes. Twelve of these 
are in the 2 front chambers CHI and CH2, where they 
are oriented along the u and v directions at +45° and 
—45° relative to the spectrometer dispersive direction 
(x); each chamber has the configuration vvvuuu. In the 
back chambers CHS and CH4 the configuration is uuvvxx 
and uuuvvv, respectively, where x indicates that the x- 
coordinate is measured, and the straws are oriented along 
the y direction. 

The individual straws are thin walled mylar tubes con- 
taining the anode wire at their center, and the gas mix- 
ture. They are made by wrapping an inner 10 /im thick 
aluminum foil, and two 50 fiin thick mylar foils around 
a mandrel, together with a heat setting glue. The inner 
diameter of the straw is 10.44 mm. The ground con- 
nection to outside is made with silver epoxy to a brass 
ferrule inserted at both ends. As shown in Fig. |3| a del- 
rin feed-through inserted in the ferrule provides gas feed 
and exhaust, and a positioning hole for a brass slit pin. 
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FIG. 3: Schematic showing the end assembly of straw tubes. 



into which the high voltage anode wire is soldered under 
prescribed stress to insure that the gravitational sagging 
and electrostatic deflection of the wire arc small. The an- 
ode wire is gold plated 25 fxm. diameter tungsten-rhenium 
wire. 

The two front chambers are identical to one another 
and contain 1008 straw tubes each. In these two cham- 
bers the straw tubes are precisely spaced by inserting 
their ends into aluminum blocks in which holes of diam- 
eter 10.75 mm have been drilled with 10.95 mm spacing 
center-to-center. Each block has 3 layers of such holes 
separated vertically by 9.5 mm and shifted by half a hole 
separation, providing a very tight packing. Each block 
accommodates 16 straws in each of the 3 layers, for a 
total of 48. The spacing between the straws of each 
plane arc maintained with mylar shims glued every 30 
cm along the length. The active area for the front cham- 
bers is 60x209 cm^. The nominal distance between the 
two front chambers is 120 cm center to center; the inter- 
vening space was occupied during this experiment by the 
100 cm gas Cerenkov detector: although not used in this 
experiment, it contributed an additional 3 mrad to the 
multiple scattering at the lowest proton energy. 

The two back chambers contain a total of 3102 straws. 
Each chamber contains six planes of straws, with the suc- 
cessive layers of straws glued together using precision 
guiding plates and pins to insure accurate positioning, 
and is protected by 0.36 mm thick carbon fiber panels at 
the top and the bottom. Both chambers arc positioned 
on a 1.9 cm thick and 31.5 cm wide plastic honeycomb, 
aluminum faced composite, that also provides a mount- 
ing surface for gas, high voltage distribution, and readout 
boards. Chambers CH3 and CH4 have active areas of 124 
x272 cm^ and 142 x295 cm^, respectively. The distance 
between these two chambers is fixed and equal to 38.0 
cm, center to center. 

The gas mixture used in the FPP chambers is 62 %/38 
% argon/ethane by volume. The straw chambers are op- 
erated at a high voltage of 1875 V. The drift velocity 
of electrons for this gas mixture and for this voltage is 
about 50 /um/ns over almost the entire volume of the 
tube. The efficiency of an individual straw tube for sin- 
gle track events is greater than 97 % after correction for 
the small gap between them. 



The analyzer consists of five sets of graphite plates with 
different thicknesses. Each set is made of two halves that 
can be moved separately on left and right. The plates 
are beveled at an angle of 45° so that the two halves 
overlap when closed. The plates have thicknesses of 1.9 
cm, 3.8 cm, 7.6 cm, 15.2 cm and 22.9 cm and they are 
separated by ~ 1.6 cm. The ability to vary the thick- 
ness of the analyzer is necessary to optimize the efficiency 
while maintaining the Coulomb multiple scattering angle 
within acceptable limits. The carbon thicknesses used in 
this experiment at different proton energies are given in 
Table I. 

The main contribution to small angle multiple 
Coulomb scattering originates in the analyzer, with ad- 
ditional contributions from the scattering in the S2 pad- 
dles, the straw tubes in the two front chambers and the 
air between them. Table I gives typical multiple scatter- 
ing angles for the relevant kinematics of this experiment. 
As Coulomb scattering is largely spin independent, in 
first order it does not affect the polarimetcr performance; 
however, multiple scattering smears out the nuclear scat- 
tering distribution, and will therefore result in a small de- 
pendence of analyzing power upon the analyzer thickness. 
Coulomb scattering results in a strong forward peak of 
protons which did not undergo nuclear scattering; these 
events are suppressed by requiring a minimum scatter- 
ing angle 2 to 3 times larger than the multiple scattering 
rms angle; here this minimum angle was fixed at 47 mrad 
(2.7°). 

TABLE I: Multiple scattering for the ten proton kinetic en- 
ergies, is the incident proton kinetic energy. Cthick. 
the analyzer thickness, ff^'^'' the root mean square Coulomb 
scattering angle, in FPP chambers (fpp), in the analyzer (C) 
and added quadratically (total). 



T 

Pine. 


Cthick. 


arms 

^fpp 




arms 
^total 


(GeV) 


(cm) 


(mrad) 


(mrad) 


(mrad) 


0.267 


7.62 


7.9 


16. 


17.8 


0.426 


22.86 


5.9 


20.8 


21.6 


0.639 


41.91 


3.8 


18.1 


18.5 


0.799 


41.91 


3.1 


16.1 


16.4 


0.959 


49.53 


2.8 


14.8 


15.1 


1.014 


49.53 


2.7 


13.1 


13.4 


1.156 


49.53 


2.2 


11.6 


11.8 


1.333 


49.53 


2.1 


10.4 


10.6 


1.599 


49.53 


1.8 


9.3 


9.5 


1.865 


49.53 


1.6 


8.2 


8.4 



To reduce the cost of electronics, the signal output 
of the individual straws is multiplexed in sets of eight. 
With multiplexing the maximum rate each tube can ac- 
cept safely is 100 kHz. One end of each straw is con- 
nected to the high-voltage distribution board, and the 
other to a readout board. The Rutgers University elec- 
tronics shop designed and built readout boards with 16 
parallel channels. The input signal from each straw (typ- 
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TABLE 11: The output pulse widths for the eight channels of 
the multiplexing circuitry on the chamber read out cards. 



Channel 


1 


2 


3 


4 


5 


6 


7 8 


Width (ns) 


25 


45 


35 


55 


85-90 


65 


100-105 75 



ically 10 mV in size) is coupled to ground through a 1500 
pF capacitor and fed into the input of an NEC1663 am- 
plifier. The amplifier output, a 100 mV positive signal, 
is fed into a LeCroy MVL407 quad comparator. This 
is a leading edge discriminator that gives a logical true 
when the input signal exceeds a supplied positive thresh- 
old voltage. The output of the comparator is then fed 
into pulse shaping circuitry. The readout board is di- 
vided into two halves, each of eight channels. The shap- 
ing circuitry for the eight channels gives a different width 
logic pulse for each of the channels. This allows all eight 
to be multiplexed together with OR chips into a single 
output channel, reducing the number of cables and chan- 
nels of TDC needed. To limit potential noise problems, 
ground planes are inserted within chamber readout card 
stacks, and differential output signals of amplitude 0.1 
V are generated. Level shifter boards located away from 
the chambers, near the TDCs, convert these signals to 
usual ECL levels for input to the TDCs. 

The output pulse widths for eight channels, generally 
adjusted to 1-2 % of the width, are given in Table 11. The 
identification of a wire within a group of 8 is obtained 
by decoding the information from pipeline TDCs which 
digitize both the leading and trailing edge times of the 
signal. The wire group is given by the multiplexed output 
carrying a signal; the actual position of the track requires 
decoding of the timing information to first identify the 
wire in the group and then calculate the drift time using 
drift velocity calibration data. Multiple tracks within the 
same subset of 8 wires cannot be decoded, unless the hits 
are separated by at least 250 ns. 



III. DATA ANALYSIS 

The data were analyzed with the standard Hall A anal- 
ysis program called ESPACE The output from ES- 
PACE includes histograms, two-dimensional plots and 
multi-dimensional arrays (ntuples), which are used in 
further analysis to obtain quantities of interest, the 
G Ep / G Mp ratio and the analyzing power Ay. 

The kinematic settings of this experiment, the cuts 
applied in ESPACE, the selection of elastic events, the 
reconstruction of tracks in the FPP chambers, and the 
cuts on FPP variables are described in subsection A. A 
description of the azimuthal event distribution and asym- 
metry after scattering in the carbon analyzer, and spin 
transport through magnetic elements of the hadron HRS, 
is given in subsection B. The methods to calculate the 
GepI Gmp ratio are described in subsection C, and sub- 



section D describes the FPP calibration. 



A. Kinematic Settings and Selection of Good 
Events 

This experiment was performed at ten different values 
of ; these values as well as other useful kinematical 
quantities are given in Table 111. 

The Hall A analysis program ESPACE calculates the 
position and angle at the focal plane, x^P,y-fP,6^P,(j)^P, 
for each event from the raw VDC data. The position, an- 
gles and momentum for proton and electron at the target, 
y, 9, (j), S, are then calculated using the HRS optics ma- 
trix; S is the relative momentum 6—^^^^, with p and Pc 
being the scattered particle's momentum and the central 
momentum of the spectrometer, respectively. 

TABLE III: Beam energies and spectrometer settings of the 
experiment. Et,eam is the beam energy, is the four mo- 
mentum transferred square, p^^, 9^^, Pp^ and 9p^ are central 
values of momentum and angle for the spectrometers detect- 
ing electrons and protons, respectively. 



Ebeam 






Sec 


PPc 




(GeV) 


(GeV^) 


(GeV/c) 


(deg) 


(GeV/c) 


(deg) 


0.934 


0.50 


0.675 


52.59 


0.756 


45.28 


0.934 


0.80 


0.516 


79.81 


0.991 


30.84 


1.821 


1.20 


1.193 


43.45 


1.268 


40.36 


3.592 


1.50 


2.815 


22.11 


1.463 


46.52 


3.592 


1.80 


2.656 


25.01 


1.649 


42.92 


4.087 


1.90 


3.093 


22.28 


1.712 


43.35 


4.087 


2.17 


2.950 


24.40 


1.872 


40.68 


4.091 


2.50 


2.774 


27.08 


2.068 


37.68 


4.091 


3.00 


2.507 


31.29 


2.357 


33.59 


4.087 


3.50 


2.241 


35.90 


2.642 


29.87 



The final optics matrix for each spectrometer was de- 
termined subsequent to this experiment using the proce- 
dure as described in Ref. 40] . Figs. ^ and show a com- 
parison between the distributions of the target variables 
y, (j), S, 9 obtained from the data and calculated using the 
Monte Carlo program MCEEP j^J] , for proton and elec- 
tron at Q2 of 3,5 GeV^, respectively. The MCEEP re- 
sults are normalized by a factor of about 0.85. This value 
seems quite reasonable, as in this experiment the trigger 
efficiency and the effect of boiling on target density were 
neither measured nor considered explicitly in the simu- 
lations. In addition, the BCM was not calibrated, as it 
would have been in the case of an absolute cross section 
measurement, for example. The agreement between the 
data and MCEEP results is good. Cuts were applied to 
all events for 9 tar (±65 mrad), (jytar (±32 mrad), ytar 
(±6.5 cm) and 5 (±5 %), to eliminate the events seen in 
the tails of Figs. Eland|31 

The experimental event rates for each and the one 
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FIG. 4: Comparison between the target variables y, 0, 5, S ob- 
tained from the data (dots) and calculated using the Monte 
Carlo program MCEEP (solid line) for protons detected in the 
right HRS. The MCEEP results are normalized by a factor 
of about 0.85, agreement between the data and the MCEEP 
results is good. 



calculated with the Monte Carlo MCEEP are given in 
the Table IV. The actual event rates are always lower 
than the MCEEP rates by about 15 %, indicating that 
there is no significant background. The large difference 
seen between MCEEP and experimental event rate at 
of 0.8 GcV^ is due to a physical aperture cutting 
acceptance for the open collimator in the electron arm. 
Table IV also includes the total number of good events, 
average current, computer dead time, and average beam 
polarization value for each value. 



TABLE IV: Experimental conditions. 





Expt. 


MCEEP 


total 


no. 


average 


dead 


average beam 




rate 


rate 


of events 


current 


time 


polarization 


(GeV") 


Hz 


Hz 






(M) 


(%) 




0.50 


1050 


1120 


2.0 X 


10^ 


4 


21 


0.560±0.030 


0.80 


250 


420 


4.6 X 


lO'^ 


10 


1 


0.544±0.006 


1.20 


1100 


1310 


1.6 X 


10^ 


24 


14 


0.497±0.012 


1.50 


420 


480 


9.2 X 


10^ 


10 


4 


0.483±0.021 


1.80 


330 


360 


1.7 X 


10^ 


13 


3 


0.611±0.020 


1.90 


1110 


1160 


5.6 X 


10^ 


63 


24 


0.391±0.002 


2.17 


830 


970 


3.8 X 


10^ 


78 


20 


0.385±0.002 


2.50 


540 


650 


5.6 X 


10^ 


74 


4 


0.390±0.008 


3.00 


300 


370 


3.2 X 


10^ 


88 


1 


0.395±0.002 


3.50 


140 


160 


2.0 X 


10^ 


94 


1 


0.384±0.007 
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FIG. 5: Comparison between the target variables y, if), 5,8 ob- 
tained from the data (dots) and calculated using the Monte 
Carlo program MCEEP (solid line) for electrons detected in 
the left HRS. The MCEEP results are normalized by a fac- 
tor of about 0.85, the agreement between the data and the 
MCEEP results is good. 



Selection of elastic ep events was accomplished by im- 
plementing a correlated cut on the missing energy 
and missing momentum Due to the kinematic con- 
straints of ep elastic scattering, no further cuts were 
needed to remove background events. The missing en- 
ergy Em is defined as: 

Em — Ebeam + Mp — {Ee + Ep) (8) 

where Ep is the scattered proton energy, and Mp is the 
mass of the proton. From conservation of momentum, 
the missing momentum, Pm, is defined as: 

^rn = \J Pm^ + -^m^ + -^m^ (9) 

Pm^ = Pe ■ sin ~pp ■ sinOp (10) 
Pmy ^ Pe ■ sin (pe - Pp ■ sin <Pp (11) 

Pm, = Pbeam ~ {Pe ' COS 6'e -|- Pp ■ COS Op) (12) 

where (pp) is the scattered electron (proton) momen- 
tum, 9e (Op) is the scattered electron (proton) Cartesian 
angle in the horizontal plane, and (pe {4>p) is the scattered 
electron (proton) Cartesian angle in the vertical plane. In 
Fig.|^ a histogram of the missing energy Em and missing 
momentum Pm is shown. There is a peak at Em — and 
a peak at about P,„ — IQMeV/c. The peak in the Pm 
histogram is not centered at zero because Pm is defined 
positive, and it has a larger width due to finite angular 
resolution in both the transverse (±2.0 mrad) and dis- 
persive (±6.0 mrad) directions. A radiative tail is seen 
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FIG. 6: Histogram of the missing energy and momentum Em 
and Pm, respectively. The peaks at Em ~ and at about Pm 
— 10 MeV/c contain the elastic ep events. The peak in Pm has 
a larger width than Em because of finite angular resolution. 



in the histogram up to about 125 MeV. The tail seen 
in the Pm histogram includes the radiative tail as well as 
events that are multiple scattered in windows. 

In Hall A, two methods of measuring the beam energy 
are now available, but at the time of this experiment 
neither was operational. The beam energy can be de- 
termined using the ep elastic kinematics, either from the 
measured scattering angles of the electron and proton or 
from Eq. ^ by forcing E,n = 0. Subsequent to this 
experiment, the beam energy was measured to a relative 
precision of about 1 x 10"'^ and the central momentum of 
the spectrometer was determined to the same precision; 
using this central momentum of the spectrometer and Eq. 

the beam energies were determined and are listed in 
Table III. Comparison to the beam energies determined 
from the scattering angles lead to the conservative con- 
clusion that the beam energy was known with a relative 
precision of 2 x 10"'^. 

Once an event was identified as an elastic ep scatter- 
ing, the next step was to search for a good track in the 
front and back FPP drift tube chambers. The analysis 
part for the FPP was incorporated in the main ESPACE 
program by the FPP group |3^ : this part of the pro- 
gram reconstructs position and angles, in the front and 
back FPP chambers, then calculates the polar and az- 
imuthal scattering angles, d and ip, respectively, and the 
position of the interaction point, Zfpp, in the carbon an- 
alyzer. Tracking in the chambers is done in the u and v 
coordinates separately. Tracking starts with identifying 
clusters of hits in chambers CHI, CH2, CHS and CH4. 



To determine a track for the front (back) chambers, there 
must be at least one hit in CHI and CH2 (CH3 and CH4) 
and at least three hits total in the front (back) chambers. 
The efficiency for an individual straw to detect a proton 
is about 97 %; the total number of hits in the front cham- 
bers is usually about 5 or 6, and it is 4 or 5 for the rear 
chambers (as the X-plane is not used). The total number 
of possible front (back) tracks is the number of clusters in 
CHI (CH3) times the number of clusters in CH2 (CH4). 
The straws have cylindrical symmetry so from a single 
hit it is impossible to distinguish if the proton passed to 
the left or right of the center wire (left/right ambiguity). 
For each possible track, a least-squares straight line fit is 
done for all possible combinations of left/right for each 
hit. Then, out of all possible tracks, the one with the 
smallest is selected as the good track. The only ex- 
ception is when the track with the smallest has only 
three hits; then if there is another possible track with 
more hits, it is selected as the good track if its is 
reasonable. 

The FPP chambers are aligned by using a software 
procedure. The carbon can be moved out to make a clear 
path between the front and back FPP chambers. Then 
the tracks in the front and back chambers are aligned to 
the tracks in the VDC so that the FPP has the same 
coordinate system as the VDC. This is done for each 
FPP chamber by adjusting the three positions: z which 
is distance from the VDC and the zero of the u and v 
axes, and adjusting the three angles of the FPP chamber 
(angle of the uz plane, angle of the vz plane and the angle 
of the uv plane). The actual alignment of the chambers 
was good, since in software the angles of the chambers 
are adjusted by less than a degree. 

A good track in the front and back FPP chambers is 
required to determine the polar angle (■(?) and azimuthal 
angle ((/?) of the scattered proton in the carbon analyzer. 
Next, we calculate the distance of closest approach be- 
tween the tracks from the front chambers and the back 
chambers, and at what distance, Zfpp, from the VDC the 
closest approach occurred. In Fig.[7| Zfpp is plotted for 
= 3.5 GeV^. The total thickness of carbon is 49.5 cm; 
it consists of 4 successive blocks which are separated by 
^ 1.6 cm, so one can see a peak for each block in Fig.[7| 
Chamber 3 is located at Zfpp of about 395 cm, as seen by 
the slight bump in Fig.E|from scattering in this chamber. 
For a selection of good events, a cut is placed on Zfpp 
depending on the thickness of carbon used for the given 
proton momentum. 

To reduce the false asymmetries a cone-test is applied 
for each event. An event passes the cone-test when the 
track that hit the back chambers at the measured polar 
angle would have hit the chambers for any possible az- 
imuthal angle. In Fig. [SJ the polar angle (i?) is plotted 
versus Zfpp with the cone-test applied. As expected, the 
range of accepted polar angles increases as the interaction 
point in the carbon is closer to the back chambers. 
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FIG. 7: Histogram of Zfpp for Q'' = 3.5 GeV with a total of 
49.5 cm of carbon consisting of 4 successive blocks of carbon 
with thickness of 22.9, 15.2, 7.6 and 3.8 cm, indicated with 
different patterns. 



B. Azimuthal Asymmetry, Spin Transport in HRS 
and Extraction of Gep/Gmp 

1. Azimuthal Asymmetry 

Proton polarimeters are based on nuclear scattering 
from an analyzer material like carbon or polyethylene; 
the proton-nucleus spin-orbit interaction results in an az- 
imuthal asymmetry in the scattering distribution which 
can be analyzed to obtain the proton polarization. The 
azimuthal angular distribution of the yield, /, of scat- 
tered protons is: 



I ^ lo il + Ay{§)pfPP-n), 



(13) 



where Iq is the unpolarized yield, P^pp is the proton po- 
larization vector at the FPP, and n is a unit vector nor- 
mal to the scattering plane defined as n = fc x fc'/ | fc x fc' | , 
with k and k' , the unit vectors in the direction of the in- 
cident and scattered proton, respectively; Ayl-d) is the 
carbon analyzing power. 

Fig. |51 shows a schematic of the polarimeter chambers 
and analyzer. It shows a non-central trajectory, where d 
is the polar scattering angle, and ip is the azimuthal scat- 
tering angle defined relative to the transverse direction 
in the polarimeter coordinate system. 

The detection probability for a proton scattered by 
the analyzer with polar angle ip and azimuthal angle ip is 
given by: 



FIG. 8: Plot of Zfpp versus the polar angle i9 with the cone 
test applied. 




FIG. 9: Schematic of the polarimeter chambers and analyzer, 
showing a non-central trajectory; i9 is the polar angle, and ip 
is the azimuthal angle from the j/-direction counterclockwise. 



2tt 



sin (f — P^PP 



(14) 

where ± refers to the sign of the beam helicity, P/^^ and 
Pn^P are transverse and normal polarization components 
in the reaction plane at the analyzer, respectively; P^^^ is 
not measured because it does not result in an asymmetry 
as seen in Eq. (|13|l . In Eq. (|14l) the helicity-independent 
polarization has been omitted because the polarization 
from two-photon exchange contributions to elastic scat- 
tering is expected to be negligible. Here e(i9, Lp) is an in- 
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strumental asymmetry that describes non-uniformities in 
detector response that might result from misaUgnments 
of the FPP chambers or from inhomogeneities in detector 
efhciency. The non-uniformity of e(i5, ip) depends upon 
the population of events on the detectors, which in turn is 
determined by the choice of kinematics. However, recog- 
nizing that the detector response is helicity independent 
and that non-uniformities are limited to a few percent, 
we may approximate Eq. (|14() by. 



/± (e(^, ± A,(P/PP sin ^ - P/PP cos ^)) 

(15) 

and obtain 



(16) 



The instrumental asymmetry can then be described 
well by the Fourier expansion 



e{'d,(p) — ao{'d) + ai{-d) cosif + a2{'d) cos2ip 

+ /3i(7?)sin^-f /32(??)sin2(^ + ... (17) 

The angular dependence of the probability distribution 
is approximated by the normalized yields. 



A^iV±r7(^) 



(18) 



where the index i refers to a bin in ip, Aip is the width 
of the bin, is the number of events in bin i, 
is the number of protons with specified helicity incident 
upon the FPP, and 77(19) is the differential efficiency of the 
analyzing reaction, defined as the ratio of the number of 
protons scattered at angle ■& with a good track in rear 
chambers over the number of incident protons with an 
acceptable track in the front chambers. It is also useful 
to define sum and difference histograms, 



E, 



(19) 
(20) 
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FIG. 10: The azimuthal asymmetry distributions for 
Q^=0.8 GeV^. The first panel shows normalized yields F"*" 
and Y~ , the second panel shows the sum distribution E, 
the third panel shows normalized yield /E and the fourth 
panel shows the difference asymmetry spectrum D. 



whose expectation values are given by. 



(A) = ^(^(P/^''sin^,-P/PPcos^,)) (21) 

(E^) = ^ (22) 

to lowest order in e. Thus, Ei measures the efficiency 
while Di is sensitive to the transverse and normal com- 
ponents of the polarization at the FPP. 

Figures [TUl and [T^ show the normalized yields and the 
sum and difference spectra for = 0.8 and 3.5 GeV^, 
respectively. Figures ITTI and [T51 show the efficiency spec- 
tra for the same settings binned with respect to polar 



angle. Note that the widths of the scattering-angle bins 
were chosen to give approximately equal statistics. 

The first panel in Figs. ^1 and E| shows helicity- 
dependent normalized yields. In an ideal polarimeter 
these distributions would exhibit reflection symmetry 
about unity, but instrumental asymmetries perturb that 
ideal pattern, especially for the larger where the 
physics asymmetry is smaller. The second panel in each 
of these figures shows instrumental asymmetries at the 
level of several percent. Figures ^2 a-nd El show these 
instrumental asymmetries divided into bins of the scat- 
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FIG. 11: The sum E = {Y+ + Y-)/2 distributions for 
Q^=0.8 GeV^ for ten ^ bins. 



tering angle z9. The smooth instrumental asymmetry seen 
in Fig. ^2 is clearly dominated by the lowest two i? bins, 
which are sensitive to chamber misalignment. This effect 
is stronger for the lowest data. The experiment was 
performed in two periods separated by about six weeks 
with the five lower point taken in the second period; 
no alignment data were taken at that time and the data 
show that the chambers must have been slightly mis- 
aligned in the intervening time. Nevertheless, the final 
results are not affected by this misalignment because it is 
helicity independent and therefore cancels when the dif- 
ference between the two helicity states is calculated. The 
efficiency spectrum for — 3.5 GeV^ contains higher 
frequency components arising primarily from large scat- 



tering angles, as shown in the bottom panel of Fig. ^1 
that arise from edge effects. These effects do not af- 
fect the final results either because they are helicity in- 
dependent also. Furthermore, omission of the large bin 
simply increases the statistical uncertainty in recoil po- 
larization by about 5 % (relative) without altering the 
reported results. 
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FIG. 12: The azimuthal asymmetry distributions for 
GeV^. The first panel shows normalized yields Y'^ 
and Y~ , the second panel shows the sum distribution E, third 
panel shows normalized yield Y^/E and the fourth panel 
shows the difference asymmetry spectrum D together with 
sinusoidal fits and values for each curve. 
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The third panel in Figs. ^] and ^] shows normal- 
ized yields corrected for instrumental asymmetry, namely 
Y^jEi. These spectra show the reflection symmetry 
expected for an ideal polarimetcr. Finally, the bottom 
panel of these figures shows the asymmetry spectra to- 
gether with sinusoidal fits in Fig. 12 that can be de- 
scribed by: 
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subsequent sections. For the present purposes it is suffi- 
cient to observe that pure sinusoids fit the difference spec- 
tra with reduced better than unity, showing that the 
instrumental asymmetries have been eliminated. There- 
fore, analysis of the difference spectrum is insensitive to 
instrumental asymmetry. 

To demonstrate helicity independence at the focal 
plane, we show the ratio of helicity plus to minus events 
at g2 of 3.5 GeV^ for variables xfP,y^P,efP, and in 
Fig. El the ratio is equal to one within statistical un- 
certainty (note the expanded y-scale) and is constant for 
each one of the four focal plane variables, thus indicating 
that there is no helicity dependence between iV+ and N~ 
events in the focal plane of the hadron spectrometer. 
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FIG. 14: The ratio iV+/iV- at of 3.5 GeV^ the ratio is 
close to one and constant for each one of the four focal plane 
variables, thus indicating that there is no helicity dependence 
of A*'^ and A'^" events in the focal plane. 



FIG. 13: The sum Ei = {Y+ 4- F")/2 distributions for 
Q^=3.5 GeV^ for ten i9 bins. 



D, = Do cos {if + S) 



(23) 



where the phase shift S would be proportional to 
(J'pGep/Gmp in an ideal spectrometer. For the final 
analysis spin precession in the spectrometer was handled 
event-by-event using the procedures to be described in 



2. Spin Transport in HRS 

The proton spin precesses as it travels from the target 
to the focal plane through the magnetic elements of the 
HRS as shown in Fig. El and therefore the polarization 
components at the target and at the FPP are different. 
The hadron HRS in Hall A consists of 3 quadrupoles and 
one dipole with shaped entrance and exit edges, as well 
as a radial field gradient. 

The primary precession angle, xsi in Fig. El is de- 
fined as the difference between the spin and the trajec- 
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tory rotation angles in the dispersive direction and can 
be expressed as: 



xe = i{t^p~i){eB + e-9fp), 



(24) 



where 7 is Ep/Mp, Hp is the proton magnetic moment, 
{6b+0—9^p) is the bending angle in the dipole for a given 
event; 6^^ is the angle at focal plane, 9 is the angle at 
the target, and in the HRS 6b — 45° for the central tra- 
jectory. As shown in Fig. 1151 in first approximation for 
the homogeneous dipole field, the transverse polarization 
component Pt is parallel to the field and does not precess; 
the longitudinal polarization component is perpendicular 
to the field and precesses with an angle xe- 




FIG. 15; Schematic drawing showing the precession by angle 
Xe of the Pe component of the polarization in the dipole of 
the HRS. 



For the HRS, the polarization vectors at the polarime- 
ter, are related to those at the target, hP, where 

P is the polarization as given by Eqs. ©,0, and Q 
through a 3-dimensional rotation matrix, (S), as follows: 



\ 



fpp 
t 

jfpp 



Snn Snt Sni 
Stn Stt Ste 
Sin Sit Su 




where h is the electron beam helicity, understood here 
to be the value of the longitudinal beam polarization 
component at the target, described with sign and mag- 
nitude. Each one of the 9 matrix elements depends 
upon a particle's trajectory in the spectrometer, and 
therefore upon the 4 reconstructed kinematical variables 
y, 4>, 5 and 9 of the recoil proton at the target. In the case 
of elastic ep scattering, there is no helicity dependent 
polarization component P„ at the target in the single- 
photon exchange approximation, however the helicity in- 
dependent component P™'^ may be induced due to the 



two-photon exchange mechanism, but this component is 
very small. The two polarization transfer components 
P/P^ and P{^^ in Eq. H14|) for each event, are then given 
by "i^: 



P/ff = SnthPt + SnihPi and 



J fpp 



StthPt + StehP( 



(25) 



With the assumption that there is only a dipole with 
an homogeneous field in front of the FPP, then Stt = 1, 
Sni — sinxe, Sti = Snt = 0, and Eq. H25|) simplifies to: 



p/PP = hPf and P/PP = hPf 



smx9- 



(26) 



However, the homogeneous dipole model is not appro- 
priate for the HRS, as it does not take into account 
the precession in the non-dispersive direction due to 
quadrupoles, fringe fields, and radial field gradient in the 
dipole. 

A better approximation is to calculate the spin matrix 
elements Sij from the bend angles in the spectrometer, 
with several less restrictive assumptions about the mag- 
netic field configuration, (i) the longitudinal component 
of the magnetic field with respect to the particle trajec- 
tory can be neglected, and (ii) the trajectory angles in 
the dipole change linearly along the path length |43| • The 
Sij elements under these assumptions are: 



Snt 
Sni 
Stt 
Sti 



sinxe sinx0 cosx'^ + 0(x^ x x^) 
sinxe C0SX4, 
COSX4, cosx'^ 



sm 



X<p, 



(27) 



where the angles 



2 sinxe 



A(l)'^ 1 



xe 

id 



X0 = 7(a*p - !)[-('/'' 

'2(1 - cosxe) 



xe 



n + W 

f cos Xe 
sinxe 
sin xe 



(/./P)(l-cosxe) 



(28) 



represent the precession in the non-dispersive direction. 
Here (j>—(j>fP is the corresponding trajectory total bending 
angle; (p'^ is the mean angle in the dipole, and A(j>'^ is half 
of the bending angle in the dipole in the non-dispersive 
plane. In first order the later two angles depend only on 
the non-dispersive target coordinate and angle: 



= (0''|y)j/ + (0''|0)0 



(29) 



The parameters, ((/)'*|?/), i4>'^\4>), (Ac/j'^ly), and (A^'^jc^) 
are the couplings to the non-dispersive coordinate and 
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angle at the target, and cannot be measured directly; 
however they can be fitted and obtained from the data 
without using any spectrometer model. An analysis in 
also shows that for the HRS, because of the rel- 
atively small angular acceptance and weak fringe field 
effects, both being sources of the longitudinal field com- 
ponent, assumption (i) above is fulfilled with good ac- 
curacy. Comparison with the full COSY calculation de- 
scribed below, shows that assumption (ii) is also a good 
approximation (see part IV C). 

The results presented here were obtained with spin ma- 
trix elements Sij calculated using a model of the HRS 
with quadrupoles, fringe fields, and radial field gradient 
in the dipole, for each tuning of the spectrometer and 
event by event with the differential-algebra-based trans- 
port code COSY 44]. For a given central momentum, 
the output of the COSY code is a table of the expansion 

coefficients (7*-''™"? gf the rotation matrix. The matrix 

'■J 

elements are calculated for each event using the x, 0, y, (f>, 
and 6 coordinates of the individual protons at the target 
with the form: 



5y = J2 C^-"'""x''e'y"'(j)''6P. (30) 

The contribution to the systematic uncertainty due to 
the model parameters will be discussed in part IV C. 



where AyP/P^ and AyP^^P are the effective asymmetries 
measured at the focal plane, and Sij are the mean values 
of the spin matrix elements; they are averages over the 
kinematical acceptance of the experiment. These linear 
equations must be inverted to obtain the target polar- 
ization asymmetries hAyPi and hAyPt, which requires 
the determinant of the averaged spin transfer coefficients, 
shown below: 



jfpp 



Snt 



Ste 



(32) 



When the precession angle xe is close to 180°, the 
determinant is close to zero as the matrix element 
Sn£ — sinxe cosx4>, and the spin transfer coefficients 
Snt and Ste are very small; hence the calculation of the 
polarization components at the target is less accurate; 
this occurs for of 1.8, 1.9, 2.5 GeV^ and more so for 
g2 of 2.17 GeV2. 

To overcome the difficulty associated with xo near 
ISO'^, a new method was developed to analyze the data; 
this method calculates the mean values of the asymme- 
tries at the target rather than at the focal plane, as de- 
scribed in Appendix A 2. Assuming, again for simplic- 
ity, that the efficiency e{ip) is constant, the determinant 
for the set of equations H51I I52fl obtained with the new 
method in Appendix A 2, is: 



C. Extraction of Ratio Gep/Gmp 

Results of this experiment published previously l45l 
were obtained using the analysis method of Ref. |43| . 
we will call this the "old method". The data presented 
here result from a complete re-analysis following a "new 
method" . In Appendix A, both methods are described 
in some detail. 

In both methods, the new and the old, the FPP data 
were analyzed in bins of polar scattering angle such 
that statistics are approximately the same in all bins. 
The analyzing power Ay appearing in formulas in the 
Appendix A and in this section is an average analyzing 
power in each d bin. The ratio Gep/Gmp was extracted 
for each bin and the ratio for a given is the average 
over all ??-bins. 

With the old method one calculates first the average 
asymmetries at the focal plane and then transports them 
to the target using the average spin matrix elements, as 
discussed in Appendix A 1. Assuming that the detector 
efficiency t{Lp) is a constant, and because the azimuthal 
scattering angle (p in the analyzer is independent of the 
target quantities x, 9, y, and 5, Eqs. (jTZI) and 
from Appendix A 1 can be written simply as: 



A P 
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(31) 
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C2 
•^tt 



o2 



SI 



SntSnl + SttSte 



(33) 



Now the determinant is non-zero even when Sne=0, 
because ^0, even ifxe — 180°. 

Both methods give very similar results, except for 
the region of xe ~ 180°, where the old method gives 
larger uncertainties, and a slightly different result but 
still within the statistical uncertainty at the of 1.8, 
1.9 and 2.5 GeV^. With the new method it became pos- 
sible to obtain a ratio value at — 2.17 GeV^, which 
has Xe — 180° for the central proton momentum. 

The values of the two asymmetries at the target, hAyPt 
and hAyPe, obtained as discussed above, can now be used 
to calculate the ratio Gep/Gmp as shown below: 



G 



Ep 



hAyPt {E^ + Ee') 



Gmp 



hAyPe 



2M 



tan — 
2 



(34) 



As can be seen from this equation, neither the beam 
polarization nor the polarimeter analyzing power need to 
be known to extract the ratio, which results in very small 
systematic uncertainties. 



D. FPP Calibration 

Elastic ep — > ep scattering provides a method to mea- 
sure the analyzing power of the analyzer and there- 
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fore calibrate the FPP. After having obtained the ratio 
Gep/Gmp in ^ given kinematics following the method 
described in section III C, the values of the polarization 
observables at the target can be calculated from Eqs. I@I 
0) rewritten as: 



Pt 



(Gep/Gmp) tan ^ 
iGlp/Gl,^ + r/e) 



Mr, 



tan^ 



(G2 



Ep 



/G\ip 



(35) 
(36) 



These are values of Pt and Pi at the target, averaged 
over the acceptance of the detectors. The asymmetry 
values hAyPt and hAyPg obtained from Eqs. 1)511 152(1 
in Appendix A 2, are then used to obtain the analyzing 
power. For each proton energy, the average value of the 
beam helicity, h from Table IV, was used to obtain the 
analyzing power Ay. 

The important properties of the FPP are the analyzing 
power Ayi^-d) and the efficiency of the analyzing reaction, 
77(19). The analyzing power ^^(1?) is plotted as filled tri- 
angles in Fig. 1161 versus the polar scattering angle at 
the ten proton energies of this experiment; for proton en- 
ergies 0.879 and 0.934 GeV the data have been placed in 
the same panel, as filled triangle and circle, respectively. 
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FIG. 16: Analyzing power as a function of polar scattering 
angle i9 in FPP, for the JLab Hall A polarimeter (filled triangle 
and circle); also shown are fits as solid line and data from 
other earlier calibration experiments as open squares. 



In each panel the energy is that of the proton corrected 
for energy loss to half the thickness of the C-analyzer 
given in Table V. For comparison with the world data, 
fits adapted to the actual thickness of the analyzer used 
in this experiment are shown. At the two lowest proton 
energies the fits are from the calibration of Waters et al. 
|32| , and in the next two they are from the calibration of 
Ransome ei aZ. 31]. All others are from the 14 parameter 
fit of the Saclay calibration of Cheung et al. [s^. In 
the same figure we also show as open squares data from 
various sources when available; in panel one the data at 
0.225 GeV are from Ref. [s^], in panels two and three the 
data at 0.440 and 0.691 GeV are from Ref. [p, and in 
panels for 1.045 GeV to 1.795 GeV, the data are from Ref. 
[33] . One must take notice that all earlier data were taken 
with significantly thinner C-analyzers. Overall the new 
data are in good agreement with fits and with previous 
data when available. 

The efficiency r/(??) for the ten values of is shown 
in Fig. 1171 For the 6 panels with the larger energies, 
the dependence is very much the same, as expected 
because the thickness of the C-analyzer was a constant 
49.53 cm (see Table I). A cut was applied to eliminate 
all data below 2.7° and above 50.0°. The lower cut was 
defined by the Coulomb scattering and the resolution of 
the drift chambers; the upper cut is determined by the 
size of the FPP back chambers. We include events with 
smaller and with larger scattering angles than in Ref. 33j 
because we used drift chambers instead of the multi wire 
proportional chambers used in the calibration of Bonin 
et al. and Cheung et al., and have therefore better posi- 
tion and angular resolution, and also we used larger rear 
chambers. The enhancement seen at small angles is in- 
dicative of Coulomb scattering events from protons with- 
out nuclear interaction and with scattering angles larger 
than 2.7°. Also as expected, the width of the forward 
peak is seen to decrease with increasing Tp. 

Combining the Ay{'d) from Fig. [Eland 77(1?) of Fig. [T7I 
gives the differential figure of merit (FOM) ri{-d)A^{d). 
The FOM values allow a quick evaluation of the number 
of "good" incident protons required to obtain a given sta- 
tistical uncertainty in the polarization components. The 
calibration data for analyzing power Ay, efficiency 77, and 
FOM are given in Table V. 

The statistical uncertainty in the results of this exper- 
iment is directly dependent on the analyzing power Ay 
and the differential efficiency r/, and can be written as: 



^pfpp 



AP/PP 



(37) 

where Ntot = + N-^ is the total number of incident 
protons on the analyzer. 
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FIG. 17: Differential efficiency as a function of proton kinetic 
energy, Tp^^^^^ , evaluated at the half thickness of the ana- 
lyzer. 



TABLE V: Averaged analyzing power, Ay, and integrated ef- 
ficiency riint. and figure of merit, FOM, of the FPP for the ten 
proton kinetic energies of the experiment, Tp^^^j^ • The range 
of integration for all three quantities is 2.7° < •& < 50°. Both 
the statistical (stat.) and systematic (sys.) uncertainties for 
FOM are given. 
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IV. RESULTS AND DISCUSSION 
A. Experimental Results 

Figure ITSl shows the results for the ratio fJ.pGEp/GMp 
from this experiment as filled circles, and from all exper- 



iments that took place after 1970 with same symbols as 
in Fig. n Our data points are shown with statistical 
uncertainties only; the systematic uncertainties are indi- 
cated by the black polygon. The data are tabulated in 
Table VI, where statistical and systematic uncertainties, 
as well as bin sizes are given for each data point. The 
results from this experiment differ from those previously 
published 45j in three ways: the points at 1.77, 1.88 and 
2.47 GeV^ have moved slightly as a result of the new 
analysis method described in subsection III C and Ap- 
pendix A, a new point was obtained at Q^=2.13 GeV^, 
and the systematic uncertainties are 2-3 times smaller. 
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FIG. 18: The ratio fipGEp/GMp from this experiment, filled 
circles, compared to world data, shown with the same symbols 
as in Fig. Q with the 2 data points of Ref. fl^ shown with 
filled squares. The absolute value of systematic errors from 
this experiment is shown by the shaded polygon. 



As seen in Fig. ^1 the electric form factor data from 
the past 30 or so years can be described by the dipole 
form factor, Gep/Gd ~ 1 considering the spread in the 
data from different experiments; however, the data from 
this experiment deviate from the dipole form factor sig- 
nificantly starting at of 1 GeV^ . Note that the results 
from this experiment are in ap par ent goo d agreement 
with the earlier data of Refs. \VA llsL Il4j which have 
much larger uncertainties; also considering the larger un- 
certainties, the SLAG data of Ref. are compatible 
with our results up to about of 2.5 GeV^. 
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TABLE VI: The ratio ^ipGEp/GMpi: statistical uncertainty 
(Icr). Asys is the systematic uncertainty from Table VII. 

2 

Q and Xe ^^e weighted average four momentum trans- 
fer squared and spin precession angle, respectively. AQ^ is 
half the acceptance. The last column Pt/Pi is the ratio of 
measured polarization components at the target, the relative 
uncertainty is the same as for UpGEp/GMp. 
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Several global analyses of the UpGep/Gmp ratio ob- 
tained by Rosenbluth separation method have been done. 
We show in Fig. three of these fits: the dotted line is 
the Bosted fit ,4^] . The solid and dashed lines are from a 
recent global fit by Arrington including all existing 
data, and a subset of the data, respectively. The sub- 
set of the data is a biased selection which was chosen to 
give the lowest possible ratio and therefore represents a 
lower limit for the ratio extracted from the cross section 
data alone. All fits tend to be dominated by the more re- 
cent SLAC data. Only the UpCsp/GMp ratios obtained 
by recoil polarization are shown here; the data from this 
paper, from ref. the early data of Milbrath et al. 

|23||. Pospischil et al. [2^, and a collection of measure- 
ments made in various Hall A experiments and analyzed 
by Gayou et al.ji^. We conclude from the comparison 
between Rosenbluth and polarization data shown in Fig. 
^|and from the further analysis done in ref. i that the 
two methods give definitively different results; the differ- 
ence cannot be bridged by either simple re-normalization 
of the Rosenbluth data or by variation of the polariza- 
tion data within the quoted statistical and systematic 
uncertainties. 

The new results are the first high precision measure- 
ments of the ratio Gep/Gmp above of 1.0 GeV^. The 
result for the ratio shows a systematic decrease as in- 
creases from 0.5 GeV^ to 3.5 GeV^, which indicates that 
Gep falls faster than Gmp- Iu the non-relativistic limit, 
this fact could be interpreted as indicating that the spa- 
tial distributions of charge and magnetization currents in 
the proton are definitely different. 



FIG. 19: The ratio UpGep/Gmp obtained in this experiment, 
and other polarization transfer experiments, including 23. . .23 . 
H^iUS' The curves are fits to the Rosenbluth separation data: 
one from Bosted |4^ shown as dotted line, and two from |4^ 
shown as solid line for all world data, and a dashed line for a 
selected subset of the world data. 



B. Nature of the Data 

To illustrate the stability of the polarization data pre- 
sented in the previous part, we show in Fig. I20l the varia- 
tion of the the polarization component ratio, Pt/ Pt at the 
target, as a function of the four target variables: y, </>, 5 
and 0, for Q^=3.5 GeV'^. In each one of the panels, the 
data are integrated over the full experimental acceptance 
for the other three variables; therefore the statistics in 
each panel are not independent. The filled circles are the 
results obtained with full COSY calculation; the open 
triangles are the values obtained not taking into account 
the spin rotation due to the quadrupoles and the field 
gradient and fringe field of the dipole. The scatter in the 
data (filled circles) for each variable is compatible with 
the statistical uncertainty on Pt/Pi, which are shown as 
the dashed lines below and above the average, shown as 
the solid line; a constant value of Pt/Pi over the accep- 
tance of the spectrometer in each one of the four tar- 
get variables demonstrates that the spin transport ma- 
trix elements are correct over the whole acceptance of 
the spectrometer, and indirectly, that the spectrometer 
model used in COSY is correct. The values of Pt/Pg are 
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FIG. 20: The ratio Pt / Pi at the target, plotted as a function of 
the 4 kinematical target variables for Q'^=3.5 GeV^. The filled 
circles are the results from COSY; the open triangles are the 
values obtained not taking into account the spin rotation due 
to the quadrupoles and the field gradient and fringe field of 
the dipole. Also shown is the range of statistical uncertainty 
corresponding to the uncertainty on fipGEp/GMp. 



The systematic decrease of the ratio fJ-pGEp/GMp 
shown in Fig. ^|and given in Table VI can be traced to 
the observed absolute values of Pt being systematically 
and increasingly smaller than the values calculated us- 
ing the dipole form factor Gd = (1 + /0-71)~^, and 
Gep/Gd and Gup/ ^-pGo- This is demonstrated in Fig. 
1211 where the values of Pt at the target obtained in this 
experiment are shown versus as solid circles. Also 
shown in this figure are the values of Pt calculated from 
Eqs. Q and lO with dipole form factors as open trian- 
gles, and open circles are P{^^ calculated from the data. 
It is seen from this figure that the additional polarization 
rotation introduced by the quadrupoles, and fringe and 
indexed fields in the dipole, is always small, and tend to 
decrease the magnitude of P}^^ from its value at the tar- 
get. The spin transport coefficient in Eq. \'6\\ Sfii is small 
but not zero; it would average to zero over a symmet- 
ric angular distribution of events in the HRS acceptance. 
However, the dependence of the cross section results 
in more events at small scattering angle than large, and 
that results in a non-zero average value of Sti, producing 
the small change of the observed transverse component 
at the target seen in Fig. |^ Noticeable is the abrupt 
bend for all three sets; this apparent discontinuity is the 
result of the choice of beam energies: the last five points 



were taken at constant beam energy, whereas the first 
five where taken with an increasing energy to optimize 
the experiment. 
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FIG. 21: The transverse component of the polarization at 
the target, Pt measured in this experiment is shown as solid 
circles. The open triangles show the value of Pt calculated 
using the dipole approximation, or scaling assumption: Ge ~ 
Gd and Gm = UpGo- Also shown as open circles are the 
values of Pt from the approximation Pt = P/^^, not taking 
into account the spin rotation due to the quadrupoles and the 
field gradient and fringe field of the dipole. 



C. Discussion of Systematic uncertainties 

The two major sources of systematic uncertainties af- 
fecting the results of this experiment arise from the align- 
ment and track reconstruction in the FPP, and from the 
spin transport in the magnetic elements of the HRS. In 
addition, from Eq. 1341 there is a small contribution due 
to the uncertainties in beam energy and scattered elec- 
tron's angle and energy; for = 3.5 GeV^, the combined 
relative uncertainty is 3 x 10^'^. 

In the FPP, the instrumental asymmetries cancel to 
first order when the difference of the azimuthal angu- 
lar distribution is calculated for the two helicities as dis- 
cussed in section HI B 1 and shown in Figs. 1101 and [T^ 
hence these asymmetries do not contribute to the system- 
atic uncertainties. No correction is necessary for dead 
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time with a polarimeter of this kind, because an event 
is defined in software by the tracks in the chambers, in- 
dependently of the experimental trigger efficiency. The 
dead-time is not helicity or polarization dependent as the 
rate of scattering are the same for the two helicities, and 
independent of the polarization of the proton. 

However, we estimate that uncertainties of the order 
of 2 mrad in the azimuthal angle, contribute to the 
systematic uncertainties. This uncertainty in 93 is due 
to the alignment accuracy of the FPP coordinate system 
relative to the reaction coordinate system at the target; 
it is based on the accuracy with which the focal plane 
chambers were surveyed, and the accuracy with which 
the FPP chambers were aligned relative to the focal plane 
chambers, as explained in part III A. 

The main contribution to the systematic uncertainties 
in the Gep/Gmp ratio comes from the calculation of the 
spin precession in the spectrometer. The COSY model 
was optimized to give the same results as ESPACE when 
reconstructing the target variables 9, 5, 0, and y from the 
measured focal plane observables 6^^ , x-^^, cj)-^Pa.nd y^^. 

The ESPACE optics matrix elements are well known 
from several optical studies of the Hall A HRSs. The un- 
certainties in the ESPACE reconstruction procedure gen- 
erate uncertainties in the COSY model, and accordingly 
in the spin transport calculations. To estimate these sys- 
tematics, one would have to vary the reconstructed tar- 
get quantities and after each variation, re-optimize the 
COSY model and recalculate the spin transport matrix. 
The COSY optimization is not a unique procedure be- 
cause one can achieve similar results by varying different 
model parameters. However the geometrical approxima- 
tion introduced in Section HIB-2, given by Eq. (|27|l . 
shows that the spin transfer coefficients depend mainly 
on the trajectory bending angles Eq. (|28|l in the disper- 
sive and non-dispersive directions and much less on any 
other target or focal plane quantities separately, indicat- 
ing that the details of the COSY parameterization are 
not important. This approximation, in which the matrix 
elements are calculated from the trajectory bend angles, 
represents the trivial part of the result. 

The ratios ^ipGEp/GMp obtained with the geometri- 
cal approximation given by Eq. H27|) . and with the full 
COSY calculation, are compared in Fig. the differ- 
ence between the results from these two methods is small 
and within the experimental uncertainty shown in Fig. 
1181 Hence, the analysis of the experimental uncertainties 
can be done by examining the geometrical approximation 
given by Eq. 1)2 7() . Accordingly, the error analysis is per- 
formed by evaluating uncertainties in the ratio Gep/Gmp 
resulting from the uncertainties in the bending angles in 
the dispersive and non-dispersive directions. Table VII 
gives the relative change in ratio at each value for the 
given change in the bending angles. 

The value of the total bend angle in the dispersive 
direction cannot be measured directly; it was obtained 
from a measurement, using the fact that P/p^=0 when 
sin xe=0. With the help of approximations given in Eq. 
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FIG. 22: The ratio UpGep/Gmp vs. Q^, solid circles, obtained 
with the COSY model, compared to results obtained using the 
approximation given by Eq. 12711 . open triangles. 



(|?7jl . P,(PP can be written as: 
pfPP = sin xe (-sin x^cos x^/iP* -t- cos X4>hPi) (38) 

The large acceptance of the HRS in xo provides a map 
of Pj(PP versus X6- In Fig- IZHl the values of the nor- 
mal component of the polarization at the FPP are shown 
versus the dispersive plane precession angle xe for of 
2.17 GeV^; also shown in this figure are the values of 
the matrix elements Sni (multiplied by the average value 
of Pe) calculated with COSY. With a 45° bend angle in 
the HRS, the precession angle is xe = 180° from Eq. 
(PHl . when the spectrometer is tuned to the central pro- 
ton momentum of 1.875 GeV/c. This procedure allows 
an accurate and independent determination of the spec- 
trometer bend angle. As demonstrated in this figure, the 
data crosses zero at 178.4°± 0.7° instead of 180.0°. The 
deviation of 1.6° in precession angle, results in a 7 mrad 
deviation of the dispersive bend angle {9b + 9 — 9fP), 
using Eq. (|24|l . A value of 7 mrad was included as an 
uncertainty in the dispersive bend angle {9b + 9 — 9^p). 

From a dedicated optical measurement described in 
Ref. [50l |. the uncertainty in the bend angle in the non- 
dispersive plane, (j) ~ 4'^^ i for the central trajectory was 
found to be 0.3 mrad. The uncertainties in the angles (f)'^ 
and A^'* given in Table VII, were calculated by varying 
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FIG. 23: The normal polarization component at the FPP as a 
function of precession angle xe- Open circles are data, dashed 
line is the fit to the data and the dots are calculated with the 
COSY model. It is seen that zero crossing occur at 178.4° 
instead of 180.0°. 



the couplings 
discussed in Ref. 
given in Table VIII. 



(A(/.''|(/.), and (A(/.'^|y), as 
. The individual uncertainties are 



TABLE VII: Dominant contributions to the systematic un- 
certainty of UpGep/Gmp, which results from the variations in 
the quantities listed in the table. The total systematic uncer- 
tainties Asys given in Table IVTl are obtained from the values 
below in quadrature. 
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^sys 


GeV^ 


7 mrad 
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0.49 


-0.004 


0.001 


0.001 


-0.004 


0.006 


0.79 


-0.009 


0.001 


0.003 


-0.004 


0.010 


1.18 


-0.017 


0.002 


0.004 


-0.003 


0.018 


1.48 


-0.025 


0.002 


0.004 


-0.005 


0.026 


1.77 


-0.034 


0.002 


0.006 


-0.002 


0.035 


1.88 


-0.032 


0.002 


0.006 


-0.005 


0.033 


2.13 


-0.032 


0.003 


0.009 


-0.006 


0.034 


2.47 


0.030 


0.003 


0.010 


-0.009 


0.033 


2.97 


0.018 


0.004 


0.009 


0.003 


0.021 


3.47 


0.009 


0.005 


0.010 


-0.001 


0.014 



TABLE VIII: Composition of the contribution in column four 
of Table VII. 
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GeV^ 


0.082 


0.038 m"^ 


0.111 
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0.49 


0.001 


0.000 


-0.001 


0.000 


0.001 


0.79 


0.003 


0.000 


-0.001 


0.001 


0.003 


1.18 


0.003 


0.001 


-0.001 


0.000 


0.004 


1.48 


0.004 


0.001 


-0.001 


0.001 


0.004 


1.77 


0.004 


0.004 


0.000 


0.001 


0.006 


1.88 


0.004 


0.004 


0.001 


0.001 


0.006 


2.13 


0.006 


0.006 


0.004 


0.001 


0.009 


2.47 


0.007 


0.003 


0.005 


0.000 


0.010 


2.97 


0.006 


0.002 


0.006 


0.000 


0.009 


3.47 


0.004 


0.003 


0.007 


0.005 


0.010 



D. Discussion of Radiative Corrections 

No radiative corrections have been applied to the re- 
sults presented in this paper. External radiative effects 
are canceled by switching the beam helicity. The internal 
correction is due to hard photon emission, two-photon 
exchange and higher-order contributions. A dedicated 
calculation of the radiative correction for single photon 
emission has been done by Afanasev et al. |5l|. Their 
calculation includes radiative corrections to asymmetries 
in elastic ep scattering for experiments in which events 
are selected on the basis of the hadronic variables only, 
i.e. the four-momentum transfer is calculated from 
{P2 — Pi)^, where p2 an pi are the final and initial proton 
four-momentum respectively, and not from the photon 
momentum, so that the integration over this photon mo- 
mentum can be performed analytically. The correction 



to the polarization observables of the proton due to single 
photon emission is of the order of a few per cent and the 
relative correction on the polarization ratio is no bigger 
than 1%. Current indications are that the contributions 
due to the other two processes mentioned above, are at 
the same percentage level. 

V. COMPARISON TO THEORETICAL 
PREDICTIONS 

The fundamental understanding of the nucleon form 
factors in terms of QCD is one of the outstanding prob- 
lems in nuclear physics. So far all theoretical models 
of the nucleon form factors are based on effective the- 
ories; they all rely on a comparison with existing data 
and their parameters are adjusted to fit the data. The 
much improved quality of the data from JLab has made 
a significant impact on theoretical models. 
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There are a number of different approaches to calculate 
nucleon form factors. Vector Meson Dominance (VMD) 
models |HllHllIllIllEilE3,IHllHi explain low to in- 
termediate behavior of the form factors. Relativistic 
constituent quark models 6Qj £1,, ,62, _63, ,645 J6j jG^, 
IHM IH^ work well in the range of this experiment. 
Other quark models that predict nucleon form factors in 
this intermediate rang e include the cloudy bag model 
[70| . di-quark model |7ll F^. and QCD sum rules [73| . 
The soliton model of Holzwarth [t^I treats nucleons as ex- 
tended objects like skyrmions and predicts nucleon form 
factors for the intermediate range. Perturbative QCD 
predicts form factor values for large Q"^ [tUtI Itt I Its I . 
In the VMD approach, the photon couples to the nucleon 
via vector mesons, whereas in QCD models the photon 
couples to the quarks directly. We discuss some of these 
calculations in more detail here and compare them with 
the new data. 



A. Vector Meson Dominance Models 

Early VMD model calculations of form factors included 
the p and its higher excited states for the isovector part, 
and the and for the isoscalar part. The number 
of mesons involved in the interaction and the coupling 
constants and masses of the mesons can be varied to fit 
the data. 

In Fig. 1^ the fits done in the 1970s by lachello et al. 
[5^ and Hohler ei al. are shown as the short dashed 
line and dashed-dot line, respectively. Their predictions 
come close to the data of this experiment, because the 
data available in the 1970s are compatible with the data 
from this experiment. 

The VMD calculation of Gari and Kriimpelmann 
incorporated a pQCD constraint to extend the range 
of the VMD description. Their fit to the data existing 
in 1985 is shown as a thin solid line in Fig. |^ A new 
fit was done by Gari and Kriimpelmann after Gep was 
re-measured with smaller error bars in Ref. 15] for 
between 1.0 and 3.0 GeV^; this fit is shown in Fig. 1241 as 
a dashed line[53|. A large difference with the earlier fit is 
seen; the new data from JLab clearly favor the older fit. 

The VMD model of Mergell et al. l55| is an expan- 
sion of the original work of Hohler et al. [S^. It includes 
the new data from Ref. 16] and a "super-convergence" 
condition to constrain the behavior of the form factors 
to the QCD predicted fall-off at large Q^. In this model, 
above = 2.0 GeV^, the parameter which indicates the 
boundary between mesonic and quark degrees of freedom 
is sensitive to pLpGEp/Gup and should be tightly con- 
strained by the new data. The result of their fit is shown 
by the dot curve in Fig. |^ A recent re-examination of 
this model [5^ using slightly modified parameters suc- 
ceeds at reproducing recent neutron form factors from 
polarization experiments, but still fails by a factor of 2 to 
reproduce the G_^/Gmp ratio from this work at Q^=3.5 
GeV^ . Lomon j57| re- visited the Gari and Kriimpelmann 
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FIG. 24: A selection of theoretical calculations with VMD 
models, compared to the data from this experiment (solid 
circles) and of Ref. p3 (empty squares). The curves are 
short dashed 52], dot-dashed|53l. thin solid and d ashed (s^l. 
dottedlHl and thick soUdlsl. 



|54l ] approach and fitted the form factor data including 
data from this experiment. This model include the p, uj 
and vector meson pole contribution and pQCD con- 
straint. Lomon also incorporated the width of the p 
meson and higher mass vector meson exchanges based 
on the Hohler-Pietarinen [s^l model. All four nucleon 
elastic form factors and the ratio PpGEp/GMp fi'om this 
experiment are fitted. The new fit for the ratio is shown 
in Fig.|31as thick solid curve. Like-wise, in 2003 lachello 
and collaborators reconsidered their 1973 work '55] ; their 
model consists of a small structure besides the pion cloud 
described by VMD. They are now able to describe the 
neutron form factors better; and still retain the original 
slope for Gep/Gmp in the Q^ range of this experiment, 
although now the ratio crosses zero at Q^ > 10 GeV^ 
|58j ; their new work includes a prediction of the nucleon 
form factors in the time- like domain 1591. 



B. Constituent Quark Models 

In the constituent quark model, the nucleon consists 
of three constituent quarks, which are thought to be va- 
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lence quarks dressed with gluons and quark-antiquark 
pairs that are much heavier than the QCD Lagrangian 
quarks. All other degrees of freedom are absorbed into 
the masses of these quarks. The early success of the non- 
relativistic constituent quark model was in describin g th e 
spectrum of baryons and mesons with correct masses |7flj. 
However, to describe the data presented here in terms of 
constituent quarks, it is necessary to include relativistic 
effects because the momentum transfers involved are up 
to 10 times larger than the constituent quark mass. 

In the earliest study of the relativistic constituent 
quark models (RCQM), Chung and Coester cal- 
culated electromagnetic nucleon form factors with 
Poincare-covariant constituent-quark models and inves- 
tigated the effect of the constituent quark masses, the 
anomalous magnetic moment of the quarks, and the con- 
finement scale parameter; the prediction is shown as a 
dotted hue in 123 

Frank et al. 63j| have calculated Gep and Gmp in the 
light-front constituent quark model and predicted that 
Gep might change sign near 5.6 GeV^; this calculation 
used the light-front nucleonic wave function of Schlumpf 
[8(tI |. The light-front dynamics can be seen as a Lorentz 
transformation to a frame boosted to the speed of light. 
Under such a transformation, the spins of the constituent 
quarks undergo Melosh rotations. These rotations, by 
mixing spin states, play an important role in the calcula- 
tion of the form factors. The results of their calculation 
are shown as the thick solid line curve in Fig. 

Several calculations with the RCQM have been moti- 
vated by the data from this experiment 0, |6^ . 
Cardarelli et al. calculated the ratio with light-front 
dynamics and investigated the effects of SU(6) symme- 
try breaking. They showed that the decrease in the ratio 
with increasing is due to the relativistic effects gener- 
ated by Melosh rotations of the constituent quark's spin 
(short dot-dashed in Fig. [SBJ. In Ref. [g^, they pointed 
out that within the framework of the RCQM with the 
light-front formalism, an effective one-body electromag- 
netic current, with a proper choice of constituent quark 
form factors, can give a reasonable description of pion 
and nucleon form factors. The results of their calcu- 
lation with two different quark form factors are shown 
as the dot-dashed and dashed curves in Fig. |55l De 
Sanctis et al. 66] have calculated the ratio Gep/Gmp 
within the hypercentral constituent quark model includ- 
ing relativistic corrections: however, the slope of their 
Gep/Gmp ratio is too small by a factor of ^2. The chi- 
ral constituent quark model based on Goldstone-boson- 
exchange dynamics was used by Boffi et al. 67] to de- 
scribe the elastic electromagnetic and weak form factors. 
They compute these form factors in a covariant frame- 
work using the point-form approach to relativistic quan- 
tum mechanics. The results of these calculations are 
shown as the short dashed curve in Fig. 

More recently Gross and Agbakpe [HI revisited the 
RCQM imposing that the constituent quarks become 
point particles as Q^ — > oo as required by QCD; using 
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FIG. 25: Comparison of constituent quark model calculations 
with the data of this experiment (solid circles ) an d from Ref. 
1^ (empty squares). The curves are: dotted [GOI] . thick solid 
[63| . short dot-dashed [fi^ . dot-dashed and dashed [fifij ]. short 
dashed [g^, and thin solid loSj. 



a covariant spectator model which allows exact handling 
of all Poincare transformations, and monopole form fac- 
tors for the constituent quarks, they obtain excellent ten 
parameter fits to all four nucleon form factors (shown as 
thin solid line in Fig. I25|l : they conclude that the recoil 
polarization data can be fitted with a spherically sym- 
metric state of 3 constituent quarks. 

For inclusiveness we mention an interpretation of the 
nucleon form factors in terms of a quark cloud and con- 
stituent quarks proposed by Friedrich and Walcher j^. 
Following on their observation that all four form factors 
show an enhancement near Q^=0.3 GeV^. Their param- 
eterization results in a very good fit to the G Ep / G mp up 
to the maximum Q^=5.6 GeV^ of the data of Ref. 48]. 



C. More Theoretical Models 

The version of the cloudy bag model (CBM) used by Lu 
et al. [73] couples a pion field to the quarks inside a bag 
such that chiral symmetry is restored. This model is an 
improvement over the MIT bag model. In CBM, quarks 
are confined within a finite spherical well of radius R of 
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about 0.8 to 1 fm. The model wave function for a proton 
is the combined wave function of the individual quarks. 
As shown in Fig. [5^1 by the dashed-dot curve, the predic- 
tion is low compared to the data from this experiment, 
but this is a simplified version of the model. The authors 
suggest that one possible future development would be 
to include tttt interactions. 

Predictions of the nucleon form factors were done by 
Kroll et al. |7lj | in a framework where the nucleon con- 
sists of quark and diquark constituents. The diquark has 
a finite size and its composite nature is taken into ac- 
count by form factors. This model describes the proton 
in terms of distribution amplitudes, and the photons and 
gluons couple to the diquark as in perturbative QCD; 
hence in the limit of — > oo the diquark picture be- 
comes the hard scattering formalism of pQCD. The pa- 
rameters of the model were determined by fits to the ep 
cross sections above = 3.64 GeV^. Two sets of pa- 
rameters gave equally good fits to the ep cross sections, 
but differed markedly in their predictions for Gep and 
the neutron form factors. The diquark-quark model is 
not expected to be applicable below — "i GeV^ . A re- 
cent prediction 81] from this model is shown in Fig.l25las 
the medium dashed curve. Another study of the di-quark 
model with connection to the generalized parton distri- 
bution is discussed in Ref . "s^l . It describes the proton in 
the double distribution formalism, as a bound state state 
of a residual quark and two quarks strongly coupled in 
the scalar and axial- vector diquark channel. The work is 
currently limited to the small region (< 0.3 GeV'^). 

In Fig.[2n|we also show the QCD sum rule based pre- 
diction of Radyushkin 73] (dot curve). This approach 
is based on the quark-hadron duality concept, which as- 
sumes that the characteristics of the free or almost free 
quarks in perturbation theory are close to the analogous 
characteristics of the hadronic spectrum integrated over 
an appropriate energy region. The validity of this model 
is in the very large region. 

In the soliton model Holzwarth [t^] applies the rela- 
tivistic corrections due to recoil and incorporates partial 
coupling to vector mesons. He uses the skyrmion as an 
extended object with one vector meson propagator and 
relativistic boost to the Brcit frame. The result is shown 
in Fig. as the dashed and solid curves, correspond- 
ing to two different strength of the w-meson coupling 
strength, g^. This model describes the ratio very well 
over the range of this experiment. 



D. pQCD Predictions 

In the p QCp approach proposed by Brodsky and 
collaborators [73 the interaction is described as a 
convolution of a hard scattering amplitude and a baryon 
distribution amplitude. Up to leading order in l/Q^j 
the magnetic form factor is proportional to a^/Q* times 
slowly varying logarithmic terms, because the momen- 
tum of the virtual photon absorbed by one quark, must 
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FIG. 26: A selection of theoretical calculations compared to 
the data from this experiment (solid circles) and from Ref. 

(empty squares). The curves are, dot-dashed j70| . short- 
dashed |8l] |. and dashed and solid |74l |. 



be shared with the two other quarks through the ex- 
change of two gluons. In lowest order the absorbed vir- 
tual photon cannot induce a quark helicity flip and con- 
sequently pQCD predicts a faster decrease of than of 
Fip with increasing Q^, by a factor Q"^ [t^- The expec- 
tation is then that Q'^F^pj F\p should become constant 
at very high . In Fig. 1271 the JLab data together with 
data of Andivahis et a/.[l^, are shown as Q'^Fipj F\p. 
The data from Ref. show flattening above of 3 
GeV^. However, the data from this experiment do not 
show yet the specific pQCD dependence. 

Recently there have been two revisions of the pQCD 
prediction for the large behavior of F^. In the first, 
Brodsky [t^] argues that the pQCD motivated behavior 
of F2 must contain an extra logarithmic term from higher 
twist contributions; the 3 free parameters a, h and c of 
the expression 



F, 



2p 



1 



ip l + (QVc)ln''(l + QVa) 



(39) 



were fitted in Ref. (77| to the data presented here aug- 
mented from the data of Ref. with the result shown 
as a solid line in Fig. In the second, Belitsky et al. 
[tsI ] reiterate the fact that the spin of a massless (or very 
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light) quark cannot be flipped by the virtual photon of 
the ep reaction. For a quark to undergo spin-flip, it must 
be in a state of non-zero angular momentum with projec- 
tion I Lz \ — As a result, the standard pQCD prediction 
for F2p (namely oc Q~^) becomes modified by a logarith- 
mic term such that: 

A 



2p 



(40) 

where A is a normalization constant; A is a cutoff con- 
stant required to suppress the infrared singularity gen- 
erated by the very soft part of the quark wave function. 
Although the constant A in the expression above is not 
determined, a fit to the data of this paper augmented by 
the data of l48l| gives A = 290 MeV, and A =0.175. The 
soft physics scale of the nucleon is determined by A; its 
size is of order of the transverse quark momentum in the 
nucleon. This fit is shown as the dashed line in Fig. 1271 
Polynomial fits to all four Sachs form factors by Kelly 
|83L 1841 indicate that the data from this experiment are 
compatible with an approach to the pQCD regime fol- 
lowing dimensional scaling. 

2.0 



1.5 - 



0.5 



0.0 




0.0 



2.0 



4.0 



6.0 



(GeV^) 



FIG. 27: Q^F2t,/Fi^ versus . The curves shown are from 
Brodsky |73 (solid line), and from Belitsky et al. (dashed 
line). The data from this experiment are shown as solid cir- 
cles, from Ref. as empty squares and from Ref. IT^ as 
empty diamonds. 

In Fig. [551 the JLab data plotted as QF2p/Fip show 
a remarkable flattening of the ratio starting at 1-2 



GeV . Inspired by the results of this experiment, Ralston 
[ssl Is^ revisited the calculation of the single-quark spin- 
flip amplitude responsible for the Pauli form factor in the 
framework of QCD and conclude that if quarks in the 
proton carry orbital angular momentum, then F2p/ Fip 
should behave like , _ . rather than the well known 



pQCD prediction of ^ (Ref. |23)- In a different ap- 
proach. Miller and Frank Wf\ have shown that imposing 
Poincare invariance leads to violation of the helicity con- 
servation rule, which results in the behavior of F2p/ Fip 
observed in the JLab data. 
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FIG. 28: The data from this experiment, from Ref. [48 
[l^ are shown as QF2p/Fip vs.Q^. 



and 



E. Generalized Parton Distributions, Lattice QCD 
and Form Factors in the Time Like Region 

The generalized parton distributions (GPDs) provide 
a framework to describe the process of emission and re- 
absorption of a quark in the non-perturbative region by 
a hadron in exclusive reactions; they are universal non- 
perturbative objects describing hard exclusive processes 
induced by photon and electrons or positrons. Pre- 
cise measurements of elastic nucleon form factors pro- 
vide stringent constraints on the parameterization of the 
GPDs. Early theoretical developments in GPDs indi- 
cated that measurements of the separated elastic form 



27 



factors of the nucleon to high may shed hght on the 
problem of nucleon spin. The first moment of the GPDs 
taken in the forward limit yields, according to the An- 
gular Momentum Sum Rule '88], a contribution to the 
nucleon spin from the quarks and gluons, including both 
the quark spin and orbital angular momentum. The t- 
dependence of the GPDs has been modeled using a factor 
corresponding to the relativistic Gaussian dependence of 
both Dirac [83 and Pauli 90] form factors of the proton. 
Extrapolation of these GPDs to i = leads to the func- 
tions entering into the Angular Momentum Sum Rule, 
and an estimate of the contribution of the valence quarks 
to the proton spin can then be obtained. This approach 
was used by Afanasev js^l who, using the Gep/Gmp data 
of this experiment in the framework of GPDs, concluded 
that valence quarks contribute about 50 % of the nucleon 
spin. When combined with inclusive deep inelastic scat- 
tering data from SMC (9^ , this result implies that about 
25 % of the proton spin comes from the orbital angular 
momentum of the valence quarks. In a recent develop- 
ment of these ideas Guidal et al have shown that the 
difhculties encountered in the Gaussian parameterization 
used in earlier work could be surmounted with a Regge 
parameterization. In one version of their parameteriza- 
tion, these authors obtain excellent fits to all 4 nucleon 
form factors. In a different approach Diehl et al use 
theoretically motivated parameterizations of the relevant 
GPDs in the very small and very large x-domains, and 
interpolated by fitting the nucleon Dirac form factors 
and F". They derive the valence contribution to Ji's sum 
rule 88]. In a related approach, Ji j|95| has shown that 
the GPDs provide a classical visualization of the quark 
orbital motion. 

There are also lattice QCD calculations predicting the 
contribution to the proton spin coming from angular mo- 
mentum of the valence quarks. For example, Mathur et 
al. jo^ have calculated the quark orbital angular mo- 
mentum of the proton from the quark energy-momentum 
tensor form factors on the lattice, and found the total 
contribution to the proton spin from the quarks to be 
60 %, of which 35 % originates from the orbital angu- 
lar momentum and 25 % from the spin of the quarks. 
Most recent calculations of the GPDs from quenched lat- 
tice QCD give u and d valence quark contributions to 
Ji's sum rule comparable to the results in Ref. [o^. A 
quenched lattice QCD calculation by Ashley et al [o^] 
with extrapolation to the quark mass gives good agree- 
ment with Gep up to Q^=l GeV^; it may be a harbinger 
of the quality of future lattice calculations. In contrast 
with all the models discussed above, lattice QCD will 
ultimately predict form factors from first principles. 

In a more general framework, the asymptotic behavior 
of both space-like and time-like elastic form factors are 
connected. Application of the Phragmen-Lindelhof theo- 
rem to the form factors shows that their asymptotic 
behavior must be the same, and that their ratio should 
go to 1 as ^ 00. The existing time-like proton form- 
factor data show that this condition is far from being sat- 



isfied at Q^=13 GeV^ as discussed by Tomasi-Gustafsson 
and Rekalo 1221 ■ However, no form factor separation has 
been possible to this point, and therefore comparison of 
asymptotic behavior in the time- and space like domains 
are currently limited by the data base. 

F. Charge and Magnetization Distribution 

Recently Kelly |83| has developed a relativistic pre- 
scription to relate the Sachs form factors to the nucleon 
charge and magnetization densities which is consistent 
with pQCD at large and with the Lorentz contrac- 
tion of the Breit frame relative to the rest frame. The 
model dependence of the fitted densities is minimized 
by using an expansion in a complete set of radial ba- 
sis functions. Details of the fitting procedure and data 
selection are in Ref. [s^]; for > 1 GeV^ the Gep 
analysis relied upon the recoil polarization data from 
Jefferson Lab and omitted Rosenbluth separation data. 
Fig. 1^ compares the fitted charge and magnetization 
densities, where both densities are normalized such that 
/o°° pi''') f^dr — 1 and where the bands include both 
statistical and incompleteness errors. The charge density 
is significantly broader than the magnetization because 
Gep is softer than Gmp, falling more rapidly with respect 
to Q2. 
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FIG. 29: The charge and magnetization distribution obtained 
by |i8S] from the data of this experiment and of Ref. [i^. 



VI. CONCLUSIONS 

We have presented the results of a measurement of 
the proton elastic form factor ratio, Gep/Gmp, oh- 
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tained in a polarization transfer experiment at JLab. 
The ratio values decrease systematically with increasing 
four- momentum transfer squared, Q^, starting at about 
1 GeV^ and up to 3.5 GeV^. Comparing the new data 
with the database of all previous experiments, we con- 
clude that the interpretation of the previous database, 
widely accepted in the literature, that the Gep/Gmp ra- 
tio is approximately constant up to 5 GeV^ , is no longer 
sustainable. As the discussion in this paper explains, 
the great strength of the polarization transfer method is 
the tight control of systematics it affords. Previous ar- 
guments about the difficulty to separate Gep from cross 
section data were based on the expectation that the ratio 
was essentially l/r/ip; but the polarization results indi- 
cate that already at 3.5 GeV^ the ratio has decreased to 
60 % of this expectation. The consequence is of course 
that at larger than 3.5 GeV^, the power of the polar- 
ization method will become even more obvious. A con- 
tinuation of this experiment to larger Q^, using the same 
technique and instrumentation, has recently been pub- 
lished jig; it shows that the ratio continues to decrease, 
reaching a value of 0.27 at Q'^='o.& GeV^ (see Fig. [TH|) . 

The re are ne w ratio results from two JLab experi- 
ments |lO0lllOl| obtained with the Roscnbluth separation 
method; these result s ag ree with previous Rosenbluth re- 
sults HI; and confirm the discrepancy be- 
tween the ratios obtained with the Rosenbluth separation 
method and the recoil polarization method. The origin 
of this discrepancy is currently the subject of intense dis- 
cussion. A likely explanation is the two-photon exchange 
process, which affects both cross section and polariza- 
tion transfer components at the level of a few percents. 
However, the two-photon effects have drastic effect on 
the results from a Rosenbluth separation method, but 
modifies the ratio obtained with polarization method by 
a few percent only. A new experiment at JLab |l02j will 
investigate the two-photon effects in the near future. 

Theoretical calculations of the proton electromagnetic 
form factors have a long history. The database for three 
out of the four nucleon form factors is limited to 10 GeV^ 
or less, reaching 31.2 GeV^ only for Gup- The basic 
physics of the interaction of the electromagnetic probe 
with the nucleon is in the difhcult region of transition 
between pictures of the nucleon, as a small core sur- 
rounded by a meson cloud, and a system of three valence 
quarks accompanied by gluons and quark-antiquark pairs 
described by QCD. At the lower end of this range, the 
assumption that the photon interacts predominantly via 
an intermediate vector meson has been very successful; 
recent reconsideration of this model provides a quantita- 
tive description of the data for all four form factors. Be- 
low of 10 GeV^, one must use non-pcrturbative QCD, 
and only QCD-based phenomenological models of the 
nucleon are available. The most successful QCD based 
model is the relativistic constituent quark model, which 
describes the drop-off in the ratio Gep/Gmp observed in 
this experiment. At a very large, but not quantitatively 
defined Q^, a perturbative version of QCD (pQCD) pio- 



neered by Brodsky and collaborators [tsI should be 
valid. An important consequence of pQCD is hadron 
helicity conservation; in terms of the non-spin flip and 
spin flip form factors (Dirac and Pauli), pQCD has gen- 
erally been understood to predict a faster decrease with 
Q"^ for F2p than Fip, by a factor of The data 

presented here clearly show that the ratio Q^F2p/Fip is 
still increasing monotonically up to 3.5 GeV^. Recently a 
careful re-examination of the pQCD prediction has led to 
the inclusion of a logarithmic factor and good agreement 
with the behavior of F2p/Fip reported in this paper. 
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VIII. APPENDIX A 

1. Old Method 

Using Eq. H14I) we construct the asymmetry f{(p) and 
the efficiency function E{(f) as: 



- (/+(^)-/-(^))/2 



2tt 



AyP/P'^smip-AyPiPPcosi 



E{v) - (/+(¥>) +r(¥'))/2 = 



27r 



(41) 
(42) 



In the next step, using Eq. (|25|l we replace P^^^ and 
Pp'^ to obtain following relations: 



27r 



hAyPtStt sini^ + hAyPiSu sin 93 
hAyPtSnt cos ip - hAyPeSni cos (p 

4^ IhAyPtXtiifi) ~ hAyPiXiiifi)] (43) 
Ztt L 



where, 
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Xt{ip) = - Snt cos If + Stt sin If 
Xe{ip) = Snecosip - S a sin ip. 



(44) 



The Fourier analyses of the asymmetry f{if) leads to the 
following expressions: 



27r 



,f{ip) cos ipdip 



hAyPt I E (if) Xt{ip) COS ipdif 



hAyPi / E{(p)Xi{(p) cos (pdip, 



(45) 



and 



2ir 



f{ip) sin Lpdip 



hAyPt 



hAyPe 



E{ip)Xt ((/?) sin ipdip 



E(ip)X{ {ip) sin Lpdip. 



(46) 

To obtain polarization components at the target, the 
integrals above are approximated by corresponding sums 
over the observed events: 



2. New Method 

In the new method we multiply the asymmetry /((/?) 
in Eq. (|43|l by Xi{ip) and Xt{pi), instead of multiplying it 



by cos(/3 and simp as in the old method, and integrate it 
over (p as shown below: 



2tt 



f{p)Xi{ip)d(p = 



= hAyPt / E{ip)Xt{(p)Xi{ip)d(p 



(49) 



hAyPe / E{ip)XUv)dip, 



and 



27r 



f{(p)Xt{(p)d(p 



hAyPt 



E{^)XU^)d^ 



hAyPe I E{ip)Xi{p)Xt{ip)dip. 



(50) 



As in the old method, the polarization components at 
the target can then be obtained by replacing the inte- 
grals in Eqs. H49I5U|) with corresponding sums over the 
observed events: 



^±cos(y3 = hAyPt At (ip) cos (p 

ev ev 

- hAyPeJ2M'P)(^osp> (47) 

ev 

zb sin (f — hAyPj ^t{^) siii(/? 

ev ev 

- AyPiJ2M^)sin'p, (48) and 



In the above event sums, the spin matrix elements Sij 
that enter in Xg^ti'p), and ip are different for each event; 
they are calculated event by event and they are functions 
of the target quantities x, 9, y, (f> and 5. 



hAyPt'^Xt{^)Xl{p) 

ev 

hAyP,Y,Xl{ip), (51) 



= hAyP^Y^Xlip,) 

ev 

- hAyPl'^Xt{p)Xt{p). 



(52) 



This matrix equation gives the target quantities hAyPi 
and hAyPt. 
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